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the  worst.  For  several  wavelengths,  a  number  of  simulative  compu¬ 
tations  are  carried  out  solving  the  equations  3TR(Ejj)/3Ad  =  0, 

3TR (Z^) /3d^^j  =  0  and  3TR(Z^)/3W  =  0  =  dial  reading  weight 

off.  Ad  =  “  ^^off^i'  ^  “  calibrating  weight)  'for  finding 

effective  distribution  of  observations  and  the  best  weights,  as 
well  as  the  worst.  A  set  of  numerical  solutions  for  the  equations 
over  a  certain  range  of  observations  is  obtained. 

Based  on  the  simulative  studies,  the  concepts  of  phase  distri¬ 
bution  and  effectiveness  of  observations  in  the  periodic  error 
calibration  are  presented  amd  so  a  design  for  the  most  effective 
distribution  of  observations  is  introduced.  For  the  calibration 
of  periodic  errors  with  several  wavelengths,  it  is  preferable 
to  select  two  weights  that  can  be  mutually  compensated  in  fitting  .. 
them  with  all  involved  periods.  Some  possible  compensative 
weights  for  an  LCR  "G"  gravity  meter  with  periodic  errors  of  1206, 
1206/17,  1206/34,  and  1  counter  units  (c.u.)  are  presented. 

'An  attempt  is  made  to  answer  how  many  observations  should  be 
made  for  determining  the  periodic  screw  errors  with  reasonable  , 
accuracy. 
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1 .  INTRODUCTION 


It  has  been  verified  by  studies  and  experiments  within  the  last  10  years 
that  the  periodic  screw  error  is  one  of  the  major  instrumental  error  sources 
of  the  Lacoste-Roraberg  (LCR)  model  "G"  and  model  "D"  gravity  meters  (Torge  and 
Kanngieser,  1979,  Kanngleser  and  Torge,  1981,  Lambert  and  Liard,  1981,  Becker, 
1981).  The  magnitudes  of  the  errors  range  from  a  few  ugal  for  short 
wavelengths  to  hundreds  of  ugal  for  long  wavelengths  (Kanngieser  and  Torge, 
1981).  It  has  become  important  to  determine  these  periodic  screw  errors  in 
laboratory  calibration,  and  then  to  eliminate  or  weaken  the  errors  while 
processing  the  data. 

An  apparatus  called  "Cloudcroft,  Jr."  is  used  in  laboratory  calibration 
(Krieg,  1981,  p.  29).  The  apparatus  Includes  two  weights.  The  larger  of  the 
two  is  referred  to  as  the  I80  mgal  weight  and  the  smaller  one  is  the  20  mgal 
weight.  Using  this  apparatus,  changes  of  approximately  20,  180,  and  200 
counter  units  can  be  induced  in  the  gravity  meter  by  placing  the  weights  in 
so-called  "welght-on"  and  "weight-off"  positions.  In  one  operation  for 
weight-off  and  weight-on,  a  pair  of  observations  is  made  at  one  location 
within  the  range  of  counter  readings.  Subsequent  observation  pairs  are  made 
at  different  locations.  Usually,  the  positions  of  the  observations  are 
distributed  with  nearly  equal  spacing  over  a  range  of  counter  readings. 

Because  of  the  nature  of  periodic  functions,  the  accuracy  of  the 
determined  parameters  of  a  periodic  screw  error  depends  not  only  on  how  many 
observations  are  made,  but  also  on  what  size  of  calibrating  weight  la  used  and 
the  distribution  of  observations.  If  the  observations  were  made  with  a 
certain  special  distribution  or  a  special  size  of  weight  (for  instance,  with  a 


distribution  of  equally-spaced  intervals  or  a  weight  equal  to  an  integral 


multiple  of  the  wavelength  of  one  periodic  screw  error),  the  errors  of  the 
adjusted  parameters  would  tend  toward  infinity.  The  accuracy  of  the 
calibration  might  be  significantly  improved  if  a  proper  distribution  is 
arranged  and  an  appropriate  size  of  weight  is  adopted.  Some  examples  will  be 
given  in  this  report. 

Determining  how  the  distribution  of  observations  should  be  designed  and 
how  the  size  of  calibrating  weights  should  be  selected  to  obtain  a  higher 
accuracy,  while  also  being  economical  of  expense  and  time,  becomes  a  problem 
of  optimization.  For  one  wavelength  of  the  periodic  screw  error,  as  a  special 
example,  a  detailed  theoretical  analysis  of  the  problem  is  presented  in  this 
report,  and  it  provides  a  basic  idea  for  solving  the  problem.  With  respect  to 
the  optimization  of  the  calibration  of  the  periodic  screw  error  involving 
several  wavelengths,  simulation  studies  should  give  ideas  of  the  most  effec¬ 
tive  distribution  of  observations  for  laboratory  calibration  (Uotila,  1982). 

A  practical  periodic  screw  error  function  for  a  LCR  model  "G"  gravity  meter  is 
given  in  the  simulation  study.  Also  in  this  report,  another  ficticious 
periodic  error  function  is  used  for  simulating  and  comparing. 

2.  MATHEMATICAL  MODELS 

2.1  Mathematical  Models  for  Adjustment 

Dial  readings  are  taken  as  the  observables.  The  function  converting 
dial  readings  to  milligal  equivalents  and  the  condition  equation  with  unknown 
parameters  are  respectively  given  by  (Uotila,  1982): 
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X  -  linear  scale  factor  (unknown  parameter) 

T.  -  pe’^iod  of  screw  error 

A,  -  amplitude  (unknown  parameter) 

> .  -  phase  lag  (unknown  parameter) 

di  -  dial  reading  (d^^  :  weight  on,  doff  '  '"'eigdt  off) 

Oj  -  mllligal  equivalent  to  dj 
W  -  calitrating  weight 

Some  other  effects,  such  as  the  earth  tidal  effect  and  drift  of  instru¬ 
ment  are  omitted  in  eq.  (2).  We  assume  that  their  effects  on  observations  are 
eliminated  in  advance.  Applying  appropriate  trigonometric  identities  (2) 
becomes 
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where  coS($j_ 

Zj^  -  Aj^  sih^j^ 


are  two  indirect  parameters.  If  A  and  p  are  directly  taken  as  parameters,  the 
linearized  equation  corresponding  to  (2)  is: 
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where  W  is  a  constant  term.  The  general  form  of  the  condition  equations  with 

parameters  is  as  follows: 

F(L  .  X  )  =  0 
a  a 

After  linearization,  the  conditions  have  the  general  form: 

BV  +  AX  +  W  =0 
o 

where  B  =  ,  A  =  =  F(L.  ,  X^) 

d  u  d  X  0  DO 

Si  Si 

The  notations  used  are  those  of  (Uotila,  1967).  From  (^)  we  get 
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The  minimum  variance  solution  of  the  system  of  linearized  condition  equations 
is  obtained  from: 

X  =  -  (a'^(bb'‘')'1a)'’a'^(bb'^)'1w 


•r- 
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.  X.  ...  B. 
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V  -  -  b’^(BB^)~''(AX  +  Wq) 


DF  -  r  -  2k  -  N 

w 

r  -  number  of  condition  equations 
Ny  -  number  of  calibrating  weights 
Concerning  the  above  expressions,  refer  to  (Uotila,  1967). 

Since  all  measurements  are  of  the  same  type,  the  weight  matrix  P  is 
assumed  to  be  the  unit  matrix.  If  x_  -  l,  then  BB^  -  2  •  1  and 

O'  p  p 

T  -1  1 

(BB  )  -  .  Therefore  we  get 
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TRd 

^g)  -  TR(ArjjA'^) 

where  -  number  of  unknowns.  Equation  (10)  can  be  used  for  checking  the 
computation.  If  Ny,  >  1  #  the  unknowns  6W  should  be  added.  Sometimes  5W  are 
neglected  because  their  effects  on  the  accuracy  of  6A  and  would  be  small. 


5 


This  has  been  tested  by  simulative  computation.  We  are  interested  only  in  the 
accuracy  of  6A  and  6(>. 


2.2  Mathematical  Models  for  Optimization 

Several  criteria  could  be  selected  for  optimization.  It  is  convenient 
for  mathematical  analysis  to  minimize  the  trace  of  that  is: 


(11)  TR(E^)  -  minimum 


Hence,  the  objective  function  for  optimization  la  as  follows: 
f  -  TR(E^)  -  f(L) 


^  “^off ^offl  •••  ^off^ 

1  2  r 


or,  -  (W,  d^^^,  Ad^ ,  Ad^,  ... 


where  Ad,  -  d  ~  d  . 

Note:  each  d^^f  is  considered  to  be  an  independent  variable.  Each  dg^  , 

J  J 

which  is  correlated  with  d^^j.  and  W,  is  a  dependent  variable. 

The  beat  distribution  of  observations  and  the  best  calibrating  weights 

are  those  which  minimize  the  objective  function.  They  can  be  determined  by 

solving  the  following  equations: 
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In  general,  -  -  TR(-5^). 
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Accordingly,  (12),  (13),  and  (14)  are 


expressed  by: 
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Based  on  differential  methods  for  matrices,  we  get 
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Then,  (15)  can  be  expressed  by; 
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It  is  easy  to  find  the  derivatives  .  There  are  three  types  of 

3L 


elements  in  matrix  A  (See  (5)),  Ad^  ,  AC .  ,  and  AS .  ^  .  Their  derivatives  with 

Oj  J  Ji 


respect  to  the  variables  of  the  objective  function  ^dj  and  dg^^ 


are  as 


Derivations  of  equations  (20)  appear  in  Appendix  A.  The  phase  lag  in  the 
above  partial  derivatives  should  be  deleted  when  -  A^cos(j)^  and  -  A^sinij)^^ 

are  taken  as  parameters.  All  of  these  formulae  for  derivatives  have  been 
checked  by  numerical  calculation  using  the  difference  method. 

Another  optimization  method,  the  so-called  "DFP"  or  "Davidon"  method,  is 
also  adopted  for  this  study.  It  is  a  modified  conjugate  gradient  method  for 
minimization.  The  steps  of  iterative  computations  are  as  follows: 
a.  select  initial  vector  for  iteration  and  a  small  positive  number  e  for 

convergence  criteria.  If  |f^  |  S  e  ,  then  the  iterative  procedure 

stops,  otherwise  proceed  to  step  b. 

-(1) 


b.  Let  k  ■  1  ,  G 


fx(X^  'M,  H 
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c.  Let  G,^ 

d.  Finding  Sj^,  put  f(X^^^)  +  into  the  minimum 

e.  Let 

f.  If  [f  ^ )  j  5  e  ,  the  procedure  stops,  otherwise  if  k  -  n  (order  of 

vector  X),  then  let  X^^^  -  and  return  to  step  a,  and  if  k  <  n,  then 

let  G^  +  i  -  fx  AX,^  -  -  X^'<\  AG,<  -  G,^+i  -  G^ 
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Let  k  ■  k  1  and  return  to  step  c.  Here  X  denotes  the  optimized  variable 

vector,  i.e.,  vector  L  mentioned  previously.  Some  constrained  conditions  are 

added  to  the  system  as  follows; 

Let  X^  ”  ( AW ,  ■SAd  ,  6Ad  ,  ...,  6Ad  ) 

1  2  r 
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or  X  -  (AW,  . 6d^^^  ) 

I  2  r 

then  the  constraints  are: 

W  S  W 

max 
^  '>»ax 

l6d^^^  I  ^  max  range  of  observations 

The  i^ethod  of  eliminated  conditions  is  used  for  processing  the  constrained 

optimization. 

Let  X,  -  Ad,  sin  Y . 

j  J  J 

W  -  4(W  ,  +  W  )  +  (W  -  W  ,  )  sin  Y  , 

n  2  min  max  max  min  n 

then  X  is  replaced  by  Y,  and  the  constrained  condition  can  be  eliminated  so 
the  problem  of  constrained  optimization  becomes  an  unconstrained  one. 


3.  OPTIMUM  DISTRIBUTION  OF  OBSERVATIONS  FOR  DETECTING  A  PERIODIC  ERROR 
FUNCTION  OF  ONE  PERIOD 


The  relative  simplicity  of  the  mathematical  and  optimization  modexs  for 
the  case  of  a  periodic  error  function  of  one  period  permits  a  rigorous 
cinalytlcal  solution.  For  error  functions  of  several  periods,  only  numerical 
analyses  are  feasible.  In  this  chapter  we  shedl  investigate  the  optimum 
distributions  of  observations  for  both  the  A,$  and  Y,Z  parameter  sets  and  the 
corresponding  optimum  weights.  The  required  number  of  measurements  to  meet  an 
accuracy  specification  for  A  or  and  certain  pathological  observation 
distributions  are  covered  as  well. 
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3.1  Optimum  Distribution  of  Observations  with  Respect  to  Parameters  Y  -  Acos^ 
and  Z  -  Asini^ 

In  this  case  the  mathematical  model  is  given  by  (3).  Considering  (6)  and 
(9)  the  trace  of  matrix  is  as  follows: 


(22) 


TR(Z^) 


2  2 
2o  X 
0 


r 

(  E 

j-1 


AS 


r 

Z 


J  J-1 


AC^  ) 


r 

Z 

J-1 


AS' 


r 

Z 

J-1 


AC 


r 

-  (  I 


sin  26 

J-1  °”J 


where 


(23)  A9 


In  accordance  with  the  discussion  in  section  2.2,  the  criterion  for 
optimization  is 

TR(Ejj)  -  minimum. 

With  regard  to  the  distribution  of  observations,  the  numerator  of  (22)  is 
constant.  Suppose  also  that  sin  AG  is  constant.  The  trace  (22)  will  then  be 
minimized  when  the  denominator  is  maximized.  This  will  occur  when  the 


conditions  hold 


0 


Z  C03  20 


Z  sin  20  -  0 

OCD 

J-1  J 


using  the  identities 


!!  _ (2j-1)ir  ,  (2J-1)Tr  „ 

Z  cos  — p; -  -  I  3in  - -  -  0,  r  *  1  , 


gives  by  comparison  with  the  above  conditions 


(24)  20 


(21  -  Dtt 


+  2mr  (n  -  1,  2,  ...). 


Substituting  according  to  (23)  and  solving  for  cl^uj j ,  the  optimum  distribution 
of  the  observations  is 

"^omj  -  I  -  1)  I  >  (J  .  1,  2 . r). 

With  the  optimum  distribution  determined,  we  seek  a  corresponding  optimum 
calibrating  weight.  Again  the  minimum  trace  condition  is  enforced.  Let  the 
distribution  related  term 

p  p 

(  Z  cos  20  f  *  {  I  sin  20  f 
J-1  J  J-1  J 

be  constant.  The  denominator  now  is  maximized  when 
2 

sin  A0  -  1 , 
from  which  follows 

A0  “  . 

Using  (23)  for  A0, 


Ad  -  d  -  d  -  4:  , 
on  off  2 

and  the  approximate  relationship 
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w  »  Ad  •  X, 

gives  the  optimum  calibrating  weight  as. 


(26)  W  -  j  X  . 


3.2  Optimum  Distribution  of  Observations  with  Respect  to  Corrections  to 
Amplitudes  and  Phase  Lags,  6A,  5i|i 

Using  now  the  mathematical  model  (U)  and  the  equations  (5)  and  (9)  the 
trace  may  be  expressed  as 

r 

a  ^  x^[(1+A^)r  +  (A^-1)  Z  cos  20„  ] 

°  1-1  1 
(27)  TR(E  )  -  - ^ - 

X  p 

A^  sin^A0(r^-((  z  cos  20  {  Z  sin  20  )^)] 

CD  ID 

j-i  J  J-1  J 

where  ~  ^ 

m.  T  offj  on^ 

The  optimum  distribution  of  measuronents  is  obtained  by  solving  the  equation 
3  TR  (Z..  ^.) 


Substituting  (27)  into  the  above  equation,  the  conditional  solutions  for  the 
equation  are  as  follows: 
a.  If  A  <  1 

A  group  of  9  whereby  the  conditions 


Z  cos  20„  -  1 

j-i  J 

(28) 

sin  20  -  0  (j  -  1 ,  2,  . . .,  r) 

"j 

pertain  will  minimize  the  trace.  From  these  conditions  it  follows 
^  ^  “n/ r+l” 

This  minimization  will  hold  when  the  number  of  measurements  is  odd.  When  r  is 


even,  the  relevant  conditions  become 
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0 


(29: 


Z  cos  20 

J-1 

sin  20 

m . 
J 


m . 
J 


(j  *  If  2f  •••#  n) 


and  it  holds 
/r-1 


r+1 


<  A  <  1 


As  A  becomes  very  much  smaller  than  the  factor/ the  criterion  (28)  becomes 
more  important  for  the  minimization  of  the  trace, 
b.  If  A  >  1 


Following  the  reasoning  above,  when  A  >  1  optimal  observations  will  be  those 
conforming  to  the  conditions 


Z  cos  20  -  -1 

J>1  “j 

(30) 

sin  20^  -  0  (j-1,2 . r) 


when  r  is  odd  and 


For  r  even,  we  require 


r 

Z  cos  20  »  0 

J-1  “j 

(31 ) 

sin  20  -  0  (j 

m  . 

J 


1.  2, 


r) 


and  therefore 


1  <  A 


to  minimize  the  trace.  The  larger  A  becomes,  the  more  negligible  the  cri¬ 
terion  (30)  is  for  minimizing  TR(Z  ,  ).  The  criteria  for  distribution  of 

6A5(^ 

observations 

sin  29  -  0  (J  -  1  ,  2,  . . . ,  r ) 

“j 

aire  necessary  conditions  for  having  minimum 

3  TR(I  )  3  TRCE,  ) 

met  simultaneously,  because  — and  -r-; - ^  have  a  common  factor 

’  '^ofr,  *  "off, 

3in29„  .  J 

With  a  derivation  similar  to  the  corresponding  one  in  section  3.1.  the 

T 

optimum  calibrating  weight  may  once  again  be  shown  equal  to  •  x. 


3.3  Requisite  Number  of  Observations  for  a  Required  Accuracy  of  A  and  4) 
The  standard  deviations  of  iX  and  6^  can  be  written  as: 


0  ^  x^(r  Z  CO320  ) 

2 _ °  J-1  _ 

6A  _  r  r 

3ln‘^A0(r‘^  -  [(  j:  00320  )  *  {  Z  3in20  )  ]} 

tn  ID 

j-1  j  j-1  j 


(32)  m 


2  2 

0  X  (r  +  Z  COS20  ) 

_ il] _ J _ 

0  P  ^  ^  2 

A  3inA0|r  -  [(  Z  CO320  )  +  (  Z  3in20  )  ]} 

ID  ID 

J-1  j  J-1  J 


Suppose  A  4C/r-1  , 

V  r*1 


and  let  i  co329 


‘  '-J.  ''u*.  *  J  •  -  *  - 


r  ■«  ' 


0 


and  z  sin20 
J-1 

Under  such  assumptions  (32)  and  (33)  reduce  to, 


(34)  m 


sinAQ  /n-1  ° 


(35)  m. 


A  slnA0  /r+1 


Alternatively,  if 


r  00320  -  I  3in20  -  0  , 

j-1  “j  j-1  ”j 


then  under  theee  conditions 


(36) 


6(() 


A  “5A 


A  3inA0  /f 


however,  13  not  exactly  minimized  in  this  case.  In  other  words, 

would  decrease  a  little,  but  would  increase  a  little  if  TR ( 5,^)  - 
minimum  is  strictly  required.  Such  regulati^^n  also  appears  in  the  optimiza¬ 
tion  for  calibration  of  more  than  one  wavelength  (see  section  4.4), 

If  the  calibrating  weight  and  required  accuracy  of  A  and  (f  are  known,  the 
number  of  observations  can  be  determined.  For  instance  assume  the  following 
data: 


Oq  -  0.002  c.u. 

T 

W  -  2  X  (x»1  mgal/c.u.) 
A  -  0 . 005  c.u. 

The  required  accuracies  aire: 
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<  0.001  tngal 

m.  <  0.087  (»  5°) 

<P 

Assuming  optimally  distributed  measurements,  to  obtain  the  required 
accuracies,  the  following  numbers  of  measurements  are  needed: 
r  >  5  for  determining  5A 
r  >  20  for  determining  5<(i 

Eq .  (35)  shows  that  as  the  amplitude  A  becomes  smaller,  the  determination 

of  the  phase  lag  <ti  becomes  more  difficult.  On  the  other  hand,  if  A  is  small 

then  the  effect  of  the  error  of  6<|i  on  the  milligal  equivalent  will  also  be 

smaller.  The  maximum  effect  on  the  error  of  iji  is  A  m  .  Letting  Am  =  0.001 

<P 

mgal ,  the  required  number  of  observations  for  both  the  determination  of  A  and 
i)  is  nearly  the  same.  For  example,  assume: 

A  *  0.005  mgal 
m^  »  Am^  a  0.001  mgal 

m^  »  _t  0.2  rad.  »  ±  11  degree 
then  r  >  5  . 

This  is  merely  a  basic  estimate  of  the  minimum  number  of  measurements  for 
the  required  accuracy  of  the  adjusted  A  and  (|>.  Below,  simulative  computations 
will  show  that  the  minimum  number  of  observations  is  approximately  45  for 

determining  a  periodic  error  function  with  7  wavelengths  if  m^^  »  Am^  <  0.001 
mgal  is  required.  However,  if  m^  is  required  to  be  less  then  5°,  the  number 
of  observations  needed  becomes  more  than  65  (see  Figs.  5a-h,  Section  4.5). 

3.4  Inadequate  Distribution  of  Observations 

In  certain  cases,  the  di st^’i bution  of  measurements  may  prohibit  the 
effective  computation  of  A  and  Several  of  such  cases  are  presented. 
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3.5  Conclusion 

The  best  distribution  of  observations  for  determining  the  unknowns  Y  -  A 

C034)  and  Z  ■  A  sin(>  is  a  uniform  distribution  within  the  range  of  one 

r  r 

period.  The  conditions  l  co329  -  I  sin2e  -  0  should  be  satisfied.  In 

j-1  j-1  ”j 

this  case,  the  highest  accuracy  for  determining  the  unknowns  A  and  would  not 
be  attained,  especially,  when  A  is  very  small.  However,  the  difference  in 
accuracy  is  not  significant. 

The  optimum  distribution  of  observations  for  the  direct  determination  of 
A  and  ^  is  also  a  uniform  distribution  over  the  range  of  several  periods.  The 
conditions  sin29^  -  0  (J  -  1,  2,  ...  r)  are  necessary  for  minimizing  both 
TR(Cjj)  ana 

The  optimum  calibrating  weight  is  [n  ♦  T  x  (n  -  1,  2,  ...) 

regardless  of  which  set  of  parameters  is  selected. 

The  TR(z:  )  and  TR(Z,  )  are  both  Inversely  proportional  to  the  number  of 
OA  6^ 

observations  r  under  the  process  of  optimization.  If  the  accuracy  for  deter¬ 
mining  A  is  required  to  be  less  than  0.001  mgad  ,  the  number  of  observations 
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has  to  be  more  than  5  when  Oq  -  0.002  mgal.  The  minimum  number  of  observa¬ 
tions  for  determining  is  20  to  attain  the  RMS  of  ij)  less  than  5*.  If  m^y^  - 
Am^^  is  required,  then  the  minimum  nunber  of  observations  for  determining  both 
A  and  is  identical. 

4.  OPTIMIZATION  FOR  DETERMINING  A  PERIODIC  ERROR  FUNCTION  WITH  MORE  THAN  ONE 
PERIOD 

4.1  Analysis  of  Factors  to  be  Optimized 

The  condition  equation  (2)  can  be  alternatively  expressed  as 

k 

(37)  Ad^x  -  2  Ad^j  3in(|i  -  ((. J  6A^  + 

k 

2  Z  Aj^  3in(|-  Ad^}  cos{^  d^  -  )  6$^^  +  W  -  0 

where,  Ad^  -  d^^  -  d^^^  -  W/x  (W  -  calibrating  weight)  d^  -  (d^ff  +  dQyj)/2. 

The  values  of  Ad^  and  sln(^  Ad^)  are  approximately  equal  to  the  size 
of  the  calibrating  weight  W  in  each  condition  equation  of  type  (37).  The 
small  value  of  the  factor  3in(^  scales  each  coefficient  of  and 

in  the  A  matrix  and  the  accuracy  of  adjusted  6Aj^  and  5(^1  will  certainly  be 
influenced  by  W.  As  is  well  known,  the  larger  W  improves  the  accuracy  of  the 
linear  scale  factor  x.  If  x  is  constant  over  the  7000  mgal  range,  then  only 
one  large  weight  is  needed  for  determining  x  with  higher  accuracy.  If  x  is 
not  constant,  an  additional  small  weight  should  be  used  for  determining  non¬ 
linear  detail  of  the  scale  factor.  In  addition,  a  larger  weight  would  improve 

the  accuracy  of  Sk^  because  there  would  be  higher  correlation  between  x  and 


The  factors  3in(|^  d  -  (bj  and  cosfl^-  d  -  (bj  are  dependent  upon  the 
distribution  of  observations.  The  manner  of  distribution  mainly  affects  the 
accuracy  of  parameters  and  6ibi  »  particularly  the  latter.  In  practice 
an  optimum  distribution  of  observations  is  referred  to  as  an  adequate  match 


between  dj^^  (or  Adj  )  and  each  period  Tj^  (i  -  1 ,  2,  . . .  k) .  Consequently  the 
general  criterion  for  optimum  distribution  is  that  observations  can  be  uni¬ 
formly  placed  over  one  period  of  each  wavelength.  The  distribution  over  the 
entire  range  of  observations  does  not  necessarily  have  to  be  uniform. 


i<.2  Optimization  of  Ccilibrating  Weights 

i<.2.1  The  optimum  weight  for  the  calibration  of  a  periodic  screw  error  with 
one  fundamental  frequency  and  second  harmonic 

As  discussed  in  Section  2.2,  optimum  calibrating  weights  are  found  by 

3  TRd  ) 

minimizing  the  objective  function  according  to  - ^ -  -  0  .  The  normal 

equations  implicit  in  are  formed  frcm  conditions  of  the  type  (37).  Zeroes 
for  the  above  minimization  problem  are  obtained  in  the  range  15  mgal  <  W  <  ^415 
mgal  and  are  determined  through  the  usual  nunerical  procedures  to  solve  a 
nonlinear  equation.  These  zeroes,  the  optimum  weights,  are  then  used  in  a 
simulative  computation  to  provide  RMS  accuracy  estimates  of  A  and  The 
following  data,  acquired  from  an  actual  bench  calibration,  was  used  for  the 
numerical  experiments: 

X  -  1  .0486J44  mgal/c.u. 

T^  -  70.94118  c.u. 

T2  -  35.47059  c.u. 

A,  -  0.0036639  mgal 
A2  ■  0.004325  mgal 


-2.08027  radians 


(p2  *  0.78761  17  radians 
r  -  1 5  observations 
doff^  -  1000.5  c.u. 

interval  -  W/x  (for  uniformly  distributed  observations). 

Two  sets  of  zeroes  corresponding  to  two  different  minima  were  found;  the 
results  are  presented  in  Tables  1  and  2.  A  set  of  weights  causing  infinite 
is  presented  in  Table  3-  Figure  1  depicts  the  behavior  of  Z^  as  a  function  of 
the  weight  W  over  the  approximate  rainge  0  <  W  <  300  mgal . 

Noting  the  ratio  WtT-jX  from  Table  1,  the  following  relationship  between 
the  optimum  calibrating  weights  and  period  T^  can  be  expressed  as 


^38)  Wopt  -  >  0.32)  T^x  (k  -  0,  1,  2,  ...) 

This  result  may  be  verified  analytically.  Suppose  that  a  weight  W  is  an  opti- 

k  IT  W  tr 

mum  weight  for  calibration  provided  it  satisfies  Z  sin  if”  7  “  Z  3in  ^  W 

i-1  i  i-1  i 

-  maximum  [W^  -  -J  .  Based  on  the  criterion,  the  formula  (38)  can  be 

proved.  Making  a  function  of  the  form:  s(W  )  -  sin  W  )  +  sin  W  }  , 

o  1,0  1 _  o 

d  S(W^) 

solve  the  equation:  — -  ”0;  if  T,  »  2T2»  we  get  a  quadratic  equation: 

0 

4  cos  ^  *  ^03  ^  -  2  =•  0  .  The  solutions  for  the  above  eq.  are 

IT  ,, 

cos  =—  W  -  - Q-  ■ 

T^  0  8 

W  . 

»  0.3183  (Taking  Positive  sign) 


- —  -  0.6817  (Taking  negative  sign) 
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ROOT  FOR  TRACE  :T,x  3in 


Table  2.  The  second  set  of  solutions  for  minimum  TR(E  ).  T,  -  70. 9*^1  c.u.  Tp  =  35. *171 


comparing  the  results  with  Table  2  and  Table  3,  the  formula  (38)  has  been 

proved,  i.e.,  Wq^  -  '''^opt*  take  the  error  of  the  solution  for 

3  TR(E  ) 

- rn -  -  0  into  account. 


4.2.2  The  optimum  weights  for  the  calibration  of  several  periodic  screw 


Based  upon  the  construction  of  the  gearbox  of  an  LCR  G  gravity  meter,  it 

may  happen  that  the  periodic  screw  error  has  more  than  two  components.  In 

this  section  the  optimum  weights  for  a  seven  component  periodic  screw  error 

3  TR(E  ) 

are  presented.  As  in  Section  4,2.1,  the  equation  - -  -  0  is  first  solved 

numerically.  The  roots  so  obtained  are  believed  to  be  in  error  by  no  more 
than  0.001  mgal.  They  are  then  used  in  simulative  computations  to  determine 
the  RMS  amplitude  and  phase  angle  data.  Results  are  presented  in  Table  4, 

In  analyzing  the  results  of  the  simulative  computations  we  may  approxi¬ 
mately  assume  that  it  is  necessary  for  an  optimum  calibrating  weight  W  to  let 
k 

E  |sin(=r-  Ad  ] I  attain  maximum.  Bearing  in  mind  that  Ad  ■  W/x, 

1-1  ‘i  °  ° 

calibrating  weights  which  maximize  this  condition  may  be  found  by  solving 


3  E  |3in(^  Ad 
J-1  ^i  ^ 

3  Ad_ 


Some  typical  results  from  such  a  computation  are  given  in  Table  5.  Using  this 

empirical  evidence  it  may  be  appropriate  to  require  the  optimum  weights 
k 

satisfy  E  j3in(^Ad  )1  >  4.3  . 


ROOT  FOR  TRACE  3in  (—  Ad  j  _ RMS 

I  MIN  Z  _  1  ° AMP 

TTl  TT"^  r^“3  r^~ii  r^“6  1  -  7  lAVE.  I  ( ugal ) 


^5  observations 


UUq 

(-W/x) 


sin  Ad^] 
i 


1  -  1 


SUM 

(3-7  only) 


18.481 

0.01 

0.05 

0.73 

1 .00 

0.88 

0.83 

1  .00 

120.502 

0.07 

0.31 

0.81 

0.95 

0.79 

0.97 

1  .00 

191.499 

0.11 

0.48 

0.81 

0.95 

0.80 

0.96 

1  .00 

23^. 'ISO 

0.14- 

0.57 

0.82 

0.94 

0.71 

1 .00 

1 .00 

262.495 

0.15 

0.63 

0.81 

0.95 

0.81 

0.95 

1.00 

26.429 

0.02 

0.07 

0.92 

0.72 

0.90 

0.80 

0.97 

44.525 

0.03 

0.12 

0.92 

0.72 

0.89 

0.80 

1  .00 

115.567 

0.07 

0.30 

0.92 

0.72 

0.87 

0.85 

0.98 

168.424 

0.10 

0.42 

0.92 

0.71 

0.91 

0.74 

0.97 

187.460 

0.11 

0.47 

0.90 

0.78 

0.63 

0.98 

0.94 

20.421 

0.01 

0.05 

0.79 

0.97 

0.96 

0.54 

0.97 

100.461 

0.06 

0.26 

0.97 

0.50 

0.72 

1  .00 

0.99 

171 .475 

0.10 

0.43 

0.97 

0.50 

0.70 

1.00 

1.00 

99.580 

0.06 

0.26 

0.95 

0.57 

0.91 

0.74 

0.97 

184.477 

0.11 

0.46 

0.95 

0.59 

0.95 

0.57 

1 .00 

Table  5.  Some  extreme  points  of  the  function  Z  sin  ^  Ad 

i-1  ^  ° 


4.2.3  Selection  of  the  optimum  weights  and  their  optimum  combination 

As  mentioned  previously,  one  large  weight  and  one  small  weight,  the  size 
of  which  would  not  be  specifically  required,  are  needed  for  determining  the 
linear  scale  factor  accurately  and  in  detail.  Likewise,  two  weights  are 


needed  for  calibrating  the  circular  error.  In  this  Instance  the  criterion 
o  ,  itW 

sin  -= -  »  sin  = —  -  maximum 


sin  -  maximum 
1 


IT  Ad, 


i3  very  important.  It  is  significant  that  the  value  sin 


has  an  effect 


on  the  accuracy  of  both  adjusted  parameters  and  4)^^.  The  ideal  condition  is 

IT  Ad 


for  the  value  of  all  sin 


terms  to  be  greater  than  0.9,  but  it  seems 


i 


impossible  to  find  one  weight  which  satisfies  the  conditions  for  all  wave¬ 
lengths.  Therefore,  we  have  to  select  more  than  one  weight  for  calibration. 
Some  of  the  optimum  weights  from  the  above  section  are  adopted  in  the  simula¬ 
tive  computation  for  further  comparison  and  optimizing. 


Numerical  experiments  demonstrate  that  the  accuracies  of  6  and  are 

IT  Ad 


very  sensitive  to  the  value  sin 


.  Table  6  shows  how  the  change  in  the 


TT  Ad, 


i 


value  of  sin 


T. 


causes  a  change  in  the  standard  deviation  of  the 


parameters  and  especially  54ij^.  If 


IT  Ad 


sin 


T: 


TT  Ad 


sin 


T. 

1 


then  m^^^  or 


is  certainly  less  than  m^^  or  m^^^  ,  respectively. 


-  Qrtj  Ul^lj 

-  A^j  »  «  A|5  =  Ay  has  been  assumed  in  the  computation  to  avoid  obscuring 

the  play  of  weight  because  the  amplitude  also  significantly  influences  the 


accuracy  of  64i<  .  Of  particular  interest  are  the  deficient  values  for 

TT  Ad 


sin 


T. 


(i  =  4,  5,  7)  for  the  180  mgal  weight,  as  evidenced  by  the  large 


RMS  values.  Similar  results  occurred  for  i  =•  3,  6  in  the  242.672  mgal 
weight . 

The  results  of  the  simulative  computation  of  the  adjustment  in  which  some 


combinations  of  weights  are  used  are  listed  in  Table  8.  They  show  that  the 

IT  Ad 


deficient  value  of  sin 


associated  with  one  weight  could  be  compensated 


1 


by  the  other  so  higher  and  homogeneous  accuracy  of  adjusted  6A^  and  6(|)^  could 


be  obtained.  It  is  obvious  that  the  value  of  RMS  6A  or  RMS  5(j)  exactly  follows 

IT  Ad 


tne  change  in  the  value  of  sin 


,  as  does  the  value  of  RMS  Y  or  RMS  Z. 


A.  r  - 

30 

observations 

wW, 

TrWg 

Period 

c.u. 

0  ' 

1 

RMS  5A 
ugal 

RMS  60 
deg 

r7 

1 

RMS  5  A 
ugal 

RMS  00 
deg 

5325 

0.10 

3.6 

0.72 

0.14 

2.8 

0.53 

1206 

0.43 

1 . 1 

6.32 

0.57 

0.7 

3.86 

70.9A1 

0.97 

0.4 

5.30 

0.73 

0.6 

7.74 

35.471 

0.48 

1  .0 

12.30 

1 .00 

0.4 

4.55 

7.882 

0.65 

1.6 

7.38 

0.90 

0.4 

5.49 

3.941 

0.99 

0.4 

5.00 

0.78 

0.5 

7.40 

1 .0 

0.89 

0.4 

5.63 

0.96 

0.4 

5.23 

Ave . 

0.64 

0.73 

PMS 

1 .5 

6.96 

1 .2 

5.45 

B.  r  - 

66 

observations 

ttWi 

RMS  641 

Period 

3ln  j— 

RMS  5  A 

RMS  60 

sin  y-r 

RMS  6  A 

c.u. 

l  .  X 

1 

Ugal 

deg 

liX 

Ugal 

deg 

5325 

0.10 

2.5 

0.51 

0.14 

1  .9 

0.34 

1206 

0.43 

0.6 

3.50 

0.57 

0.4 

2.60 

70.941 

0.97 

0.3 

3.33 

0.73 

0.3 

4.63 

35.471 

0.  48 

0.5 

6.99 

1.00 

0.3 

3.27 

7.882 

0.65 

0.4 

5.15 

0.90 

0.3 

3.65 

3.941 

0.99 

0.3 

3.15 

0.78 

0.3 

4.36 

1.0 

0.89 

0.3 

3.80 

0.96 

0.3 

3.44 

Ave . 

0.64 

0.73 

RM3 

1  .0 

4.30 

0.8 

3.45 

C.  r  - 

90 

observations 

itW^ 

RMS  6  A 

RMS  64 

Period 

3in  =— 

RMS  5  A 

RMS  60 

sin  =-;j- 

c.u. 

T^x 

ugal 

deg 

ugal 

deg 

5325 

0.10 

2.1 

0.42 

0.14 

1.6 

0.30 

1206 

0.43 

0.5 

2.99 

0.57 

0.4 

2.22 

70.941 

0.97 

0.2 

3.02 

0.73 

0.3 

3.36 

35.471 

0.  48 

0.5 

6.02 

1 .00 

0.2 

2.76 

7.882 

0.65 

0.3 

4.39 

0.90 

0.2 

3.13 

3.941 

0.99 

0.2 

2.90 

0.78 

0.3 

3.63 

1 .0 

0.89 

0.2 

3.22 

0.96 

0.2 

2.96 

AVE. 

9.64 

0.73 

RMS 

0.8 

3.64 

0.7 

2. 51 

W,  •  '80  ngai  ,  Wj  -  2‘-i2.672  mgal 
A,  -  0.2957  mgal,  A2  -  0.0095  mgal 
Aj  •  A,j  •  A^  •  A^  •  Ay  •  0.00^3  ngal 

Table  6.  TMe  Relaiibnshlp  Between  the  Accuracies  of  Adjusted  Parameters  ana 

T  W 

'/alue  bf  sin  yT-jj-  According  to  the  '/ariant  Number  of  Measurements 


■»  I  !_■  ,1  ^im  11 1.»  I  L«  ■  ■  J>T» 


I  IP*  giu^U"U"U"i."  i/'  V  ’ 


A .  W 1  * 

242.672  mgal , 

W,  -  1 30.0  mgal 

r,  -  50  obs. 

'’1  ■ 

:?0  obs . 

rp  -  1 6 

•■'2  * 

40 

ttW, 

itW^ 

Period 

sin  =—7 

sin  f—T 

RMS  6  A 

RMS  S<P 

RMS  5  A 

RMS  (5o 

c.u. 

T^x 

i  .  X 

1 

ugal 

deg 

ugal 

deg 

5325 

0.14 

0.10 

2.0 

0.38 

1  .8 

0.34 

1206 

0,57 

0.43 

0.5 

2.78 

0.4 

2 . 48 

70. 9^*1 

0.73 

0.97 

0.4 

4.94 

0,3 

3.50 

35.1^71 

1,00 

0.48 

0.3 

3.61 

0.3 

3.43 

7.382 

0,90 

0.65 

0.3 

3.91 

0.3 

3.54 

3.9*11 

0.78 

0.99 

0.4 

4.59 

0.3 

3.  41 

1.0 

0,96 

0.39 

0.3 

3.64 

0.2 

3.'2 

Ave . 

0.73 

0.  on 

RMS 

0.3 

3 .68 

0.7 

3.03 

B.  W.  - 

242.672  mgal , 

'Wp  -  38.540 

mgal 

1 

r,  -  50 

obs . 

I'l 

.  50 

obs . 

r2  •  1 6 

’'2 

-  40 

irW^ 

nW^ 

Period 

sin  =— 

Sin  =-r 

RMS  SA 

RMS  5()i 

RMS  5A 

RMS  5o 

c.u. 

T^x 

T^x 

Ugal 

deg 

Ugal 

deg 

5325 

0,14 

0,05 

2.1 

0,39 

2.0 

0.3s 

1206 

0,57 

0.22 

0,5 

2.91 

0.5 

2,31 

70.9*11 

0.73 

0.56 

0,4 

4.77 

0.3 

4.28 

35. *171 

1 ,00 

0.93 

0.3 

3,28 

0.2 

2.82 

7.382 

0,90 

0.79 

0.3 

3.73 

0.3 

3.28 

3.9*11 

0,78 

0.97 

0.3 

4.25 

0.3 

3.50 

1  .0 

0,96 

0.98 

0.3 

3,66 

0.2 

3.18 

AVE. 

0,73 

0.64 

RMS 

0.3 

3.54 

0.8 

3.10 

c.  w,  . 

.  242.672  mgal 

.  Wp  *  ’5*393 

mgal 

r^  •  50 

obs . 

f’l 

«  50 

obs . 

r2  -  16 

'  2 

-  aO 

ffW  ^ 

wWp 

Period 

sin  - — 

sin 

RMS  iA 

RMS  50 

RMS  5  A 

RMS  64 

c.u. 

T.  X 
i 

ugal 

deg 

Ugal 

deg 

5325 

0.  1  4 

0.01 

2.1 

0.40 

2.1 

0.40 

1206 

0.57 

0.05 

0.5 

2.97 

0.5 

2 . 96 

70.9*11 

0.73 

0.1'3 

0.3 

4.56 

0.3 

3.38 

35  471 

1 ,00 

1  .00 

0.3 

3.38 

0.2 

2.90 

7.332 

0.90 

0.88 

0.3 

3.65 

0.2 

3.13 

3.9*11 

0.78 

0.82 

0.3 

4.16 

0.3 

3.50 

1 .0 

0.96 

1.00 

0.3 

3.32 

0.3 

2. '6 

AVE. 

0.73 

0.64 

RMS 

0.9 

3  .  44 

0.8 

2.98 

A,  -  0. 

2957  mgal ,  A2 

-  0.0095  ingal 

A3  .  Ay 

-  A5  .  Ag  -  A 

-  0.  DO43  mgal 

Table  7.  The  Results  of  Si.-nulati '/e  Adjustment  for  Estimating  and  Eomparing 
the  Accuracies  of  Adjusted  Parameters  According  to  the  Variant 
Eorabination  of  Twc  Weights  in  the  Eaiitrat.on. 
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A.  W,  . 

130.0  mgal , 

W2  -  53.540 

mgal 

r  •  *50 

obs . 

r.  -  50  0  bs . 

r2  *16 

r2  .  4C 

ttW  ^ 

nWo 

Period 

1 

Sir. 

RMS  5  A 

RMS  5o 

RMS  6  A  RMS  60 

o.u. 

■4  gal 

deg 

ugal  ieg 

5325 

0.10 

0.05 

2.7 

0.58 

2.5  7.55 

1206 

0.43 

0.22 

0.6 

3.76 

0.6  3.54 

70.941 

0.97 

0.56 

0.3 

3.95 

0.3  3.63 

35.471 

0.48 

0.93 

0.4 

5.38 

0.3  3.52 

7.382 

0.65 

o.rg 

0.4 

4.39 

0.3  5.96 

3.941 

0.99 

0.9" 

0.3 

3.63 

0.2  3. '6 

1 .0 

0.39 

0.98 

0.3 

3.74 

0.2  3.2" 

AVE. 

0.64 

0.64 

*  «  T 

i.r 

1.0  3.33 

3.  W,  . 

180.0  mgal , 

Wp  -  20.0  mgal 

r,  -  50 

obs . 

r.  -  50  CDS. 

r2  •  1 6 

r2  -  40 

irW  ^ 

FT 

Period 

sin 

RMS  «A 

RMS  60 

RMS  6A  RMS  6® 

c.u. 

TiX 

i 

ligal 

deg 

ugal  leg 

5325 

0.10 

0.01 

2.7 

0.60 

2.7  0.60 

1206 

0.43 

0.05 

0.7 

3.90 

0.7  3.88 

70.941 

0.97 

0.75 

0.3 

3.81 

0.2  3.36 

35.471 

0.48 

0.99 

0.4 

5.37 

0.3  3.76 

7.882 

0.65 

0.97 

0.3 

4.61 

0.3  3.52 

3.941 

0.99 

0.48 

0.3 

3.77 

0.3  3.55 

1.0 

0.39 

0.22 

0.3 

4.27 

0.3  “.'9 

AVE. 

0.64 

0.50 

RMS 

1  . 1 

4.01 

1.1  3 . 45 

c.  w,  . 

1 80 . 0  mgal , 

W2  -  19.393 

mgal 

r,  -  50 

obs . 

-  50  003. 

r2  •  1  6 

r,  -  40 

wWg 

Period 

sin  ^ — 

sin  yT’ 

RMS  «A 

RMS  60 

RMS  6  A  R.MS  6$ 

c.u. 

A  .  X 

1 

i 

ugal 

deg 

ugal  deg 

5325 

O.'O 

0.01 

2.7 

0.60 

2.7  0.6C 

1206 

0.43 

0.05 

0.7 

3.91 

0.7  3.88 

70.941 

0.97 

0.73 

0.3 

3.37 

0.2  3.38 

35.471 

0.48 

1  .00 

0.4 

5.83 

0.3  4.06 

7.882 

0.65 

0.88 

0.3 

4.71 

0.3  3.72 

3.94' 

0.99 

0.82 

0.3 

3.61 

0.2  3.11 

1  .0 

0.39 

1 .00 

0.3 

3.76 

0.3  2.37 

AVE.  0.64  0.64 

RMS  ' • ’  4.04  ' . 1  3.28 

A.  -  0.2957  Jigai,  Aj  -  0.0095  rngai 
Aj  •  Au  •  A^  -  Ag  -  A^  •  0.0043  nigai 

T.iDle  S.  The  Results  of  Simulative  Adjustment  for  Estimating  and  Comparir 
the  Accuracies  of  Adjusted  Parameters  According  to  the  Variant 
Oodibi  nati  on  of  Two  Weights  in  the  Calibration. 


Concerning  the  number  of  observations,  it  is  reasonable  to  distribute 


more  observations  for  using  the  weight  which  has  the  larger  value  of  the  sum 
IT  Ad 


of  sin 


(see  Tables  7,  8).  One  large  weight  is  necessary  for  improving 


the  accuracy  of  the  parameters  of  the  long  periodic  screw  error,  and  the 


weight  could  also  be  used  in  precisely  determining  short  periodic  screw  errors 

TT  Ad 

provided  it  has  a  large  value  for  sin  — -  for  each  short  wavelength.  If 

non-linear  scale  factor  detail  is  not  taken  into  account,  it  would  be  better 

to  select  two  large  weights  which  could  be  compensated  with  each  other  in  the 
IT  Ad 

value  of  sin  -  . 

^  i  IT  Ad 

The  sufficient  compensation  of  the  defective  value  of  sin  -= -  with 

^i 

the  other  should  be  considered  as  a  principle  for  selecting  an  optimum  com¬ 
bination  of  calibrating  weights.  We  suggest  that  the  combination  of  weights 
should  be  selected  by  the  following  criteria: 


7  ttW 

(39)  I  1  sin  ~  I /7  ^  0.88 

i  ”  3  i 

ttW  iiW 

(40)  lain  ^  I  -  (sin  ^  j  >  i.6 

(i  •  3.  ...  7) 


where  W,  is  a  large  weight  and  W2  a  compensative  weight.  According  to  these 
criteria  and  Table  5,  several  possible  combinations  of  weights  are  selected  in 
Table  9.  Note  that  the  weight  of  option  depends  on  the  scale  factor  x, 
because  W  »  Ad^  x.  There  are  some  differences  in  the  scale  factor  x  between 
different  gravity  meters;  therefore,  no  universal  weight  can  be  fitted  for  all 
models.  If  a  change  in  the  value  of  x  is  assumed  to  be  0.005,  then  the 


corresponding  variation  of  weight  W  would  be  from  about  0.1  to  1.0  mgal .  The 


I  AVERA 


oorreaponding  change  in  the  value  of  ain  -= — -  can  range  from  0.15/Tj^  to 
1.5/Tj^.  It  ia  obvloua  that  the  change  will  have  a  aignif Leant  effect  on  the 
accuracy  of  determining  5A  and 

We  auggeat  that  every  gravity  meter  ahould  have  ita  own  apecial  calibrat¬ 
ing  weighta  which  fit  with  ita  scale  factor  and  perioda  of  circular  errora. 

The  apecial  calibrating  weighta  ahould  be  conaidered  aa  an  acceaaory  and 
aupplied  with  gravity  metera  to  the  uaera. 

4.3  The  Optimum  Intervala  for  Uniform  Diatribution  of  Obaervatlona 
4.3.1  Four  illuatrative  numerical  examplea 

It  ia  very  important  for  a  uniform  diatribution  to  have  an  adequate 
interval  of  obaervationa  becauae  of  the  periodicity  of  the  circular  error 
function.  We  ahould  carefully  arrange  the  interval  for  making  uniformly 
diatributed  meaaurementa.  Otherwiae,  the  reaulta  would  be  incorrect,  regard- 
leaa  of  how  many  obaervationa  are  made.  In  thia  aection  the  reaulta  of  four 
numerical  examples  are  preaented  in  the  following  pages  to  illustrate  some 
important  points.  The  observations  which  were  used  for  the  experiments  may  be 


found  in  Appendix  C. 


Experiment  1 


r  -  number  of  obeervationa  -  30 

Oq  -  a  priori  variance  of  observations  -  0.0018  c.u. 

W  -  calibrating  weight  -  180.0  mgal 
Ad  »  180  counter  units 

Adjusted  Parameters 

Scale:  x  -  1.0^8  ragal/c.u. 

A  posteriori  standard  deviation  of  scale:  3^^  -  0.0005  mgal/c.u. 


i 

Amplitude 

Phase 

Y  -  A  cos  (ji 

Oy 

Z  -  A  sin  4) 

( ugal ) 

(deg) 

(ugal ) 

(ygal) 

1 

246.83 

22.4 

228.3 

96.4 

93.9 

40.1 

2 

9.00 

238.3 

-4.7 

1 .5 

-7.7 

3.7 

3 

2.92 

254.3 

-0.8 

0.6 

-2.8 

1  .2 

4 

4.96 

70.3 

1.7 

1.5 

4.7 

1.3 

5 

1.33 

342.1 

1.3 

1.4 

-0.4 

0.9 

6 

2,60 

288.0 

0.8 

0.4 

-2.5 

1.2 

7 

44.72 

193.5 

-43.5 

60.5 

-10.4 

33.3 

Matrix  of  correlations  for 

Amplitudes 

(standard  deviations  along 

diagonal  in 

ugal) 

1  104.  0.95 

0.96 

0.82 

0.85 

0.96 

0.38 

2  3.79 

0.91 

0.73 

0.71 

0.92 

0.26 

3 

1.32 

0.80 

0.82 

0.93 

0.38 

4 

1 .69 

0.72 

0.78 

0.66 

5  symmetric 

1.39 

0.80 

0.33 

6 

1.19 

0.36 

7 

28.6 

RMS  steindard  deviation  -  Amplitude:  41 

ugal 

Matrix  of  correlations  for 

Phase  angles 

(standard  deviations  along 

diagonal  in 

degrees) 

1  0.70  -0.24 

-0.93 

0.38 

-0.28 

0.39 

-0.29 

2  6.88 

0.81 

-0.23 

0.55 

-0.82 

0.19 

3 

5.72 

-0.20 

0.21 

-0.11 

0.32 

4 

1  3.27 

-0.45 

0.1  1 

-0.54 

5  symmetric 

22.37 

-0.12 

0.49 

6 

5.70 

-0.25 

7  809.3^ 

RMS  Standard  deviation  -  Phase:  310° 
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Experiment  2 


r  -  30 

Oq  -  0.0018  mgal 
W  -  180.0  mgal 
Ad  -  1  71 . 67  G  .u . 

Adjusted  Parameters 


Scale  X  -  1.048  mgal/c.u. 
0^  -  0.00316  mgal/c.u. 


1 

Amplitude 

Phase 

Y  -  A  cos  (p 

A 

Oy 

Z  -  A  si 

(ugal) 

(deg) 

(ugal) 

(ugal) 

1 

295.72 

21 .7 

276.6 

3.7 

110.2 

2 

9.50 

242.6 

-4.4 

1.1 

-8.5 

3 

3.66 

246.2 

-1.5 

0.4 

-3.5 

14 

4.32 

41 .8 

3.2 

1  .0 

2.9 

5 

2.32 

1.4 

2.6 

0.6 

0.1 

6 

3.56 

300.6 

1 .5 

0.4 

-2.5 

7 

4.45 

257.5 

-0.9 

0.4 

-3.9 

Matrix  of  correlations 

-  Amplitude 

(standard  deviations  along  diagonal 

in  ugal) 

1  3.68  -0.02 

0.01 

-0.02 

0.01 

-0.01 

2  1.09 

-0.09 

0.50 

-0.19 

0.08 

3 

0.41 

-0.05 

-0.05 

-0.20 

4 

0.96 

-0.15 

0.04 

5 

0.60 

-0.03 

6 

7 

RMS  standard  deviation 

-  Amplitude: 

1 .5  ugal 

0.39 

Matrix  of  correlations 

-  Phase 

(standard  of 

deviations 

along  diagonal 

in  degrees) 

1  0.75 

-0.02 

0.04 

-0.06 

0.10 

0.05 

2 

6.67 

0.02 

-0.49 

0.19 

0.04 

3 

6.19 

-0.07 

0.05 

0.19 

4 

13.3 

-0.24 

-0.09 

5 

15.3 

0.06 

6 

7 

RMS  standard 

deviation 

-  Phase;  9® 

6.38 
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Experiment  3 


r  -  30 

W  -  203.89^4  mgal 
Ad  -  1 9^. ^36  G .u. 


Standard  deviations  of  adjusted  parameters 
RMS  5x  -  0.003  pgal/c.u. 


RMS 

^2 

"3 

Ti, 

"5 

^6 

h 

Mean 

RMS 

6 A  (ygal) 

3.33 

0.77 

0.54 

0.39 

0.60 

0.55 

0.41 

0.34 

6(1)  (deg) 

0.62 

4.85 

6.79 

5.21 

72.57 

79.56 

5.56 

40.93 

Y  (ygal) 

3.29 

0.78 

0.51 

0.39 

6.03 

6.01 

0.41 

3.48 

Z  (pgal) 

3.23 

0.79 

0.54 

0.39 

5.45 

5.45 

0.41 

3.19 

Experiment 


r,  -  15,  r2  -  15 

-  203.894  mgal,  W2  -  46.691  mgal 
Ad^  -  194.436  G.U.,  Ad2  -  44.525  c.u. 

Standard  deviations  of  adjusted  parameters 


RMS  5x  -  0 

.008  ugal/c 

.u. 

RMS 

Tl 

^2 

^3 

T4 

6A  (ugal) 

9.67 

1.13 

0.55 

0.41 

(deg) 

0.85 

6.62 

6.96 

7.50 

Y  (ugal) 

8.95 

1.11 

0.52 

0.41 

Z  (ugal) 

5.72 

1.12 

0.55 

0.56 

Let  U3  now  analyze  the  results  of  the  four  experiments.  In  Section  3.^4 
we  saw  that  the  standard  deviation  of  adjusted  parameters  will  tend  to 
Infinity  If  Ad  -  nT^^  (n  being  a  positive  Integer),  l.e.  the  Interval  of  obser¬ 
vations  is  equal  to  an  integral  multiple  of  period  Tj^.  In  experiment  1,  the 
standard  deviation  of  skj  -  28.6  pgal  and  the  standard  deviation  of  6<pj  -  309 
degrees  are  large  because  Ad  -  180  Ty.  In  experiment  2,  Ad  -  W/x  -  171.67  is 
not  equal  to  any  Integral  multiple  of  any  period  so  the  accuracy  is 
immediately  Improved. 

There  are  many  Inadequate  intervals  Ad^  which  are  approximately  equal  to 
nTj^.  In  this  case  it  is  easy  to  understand  the  reason  why  the  intervals  Ad  - 
nTj^  are  inadequate.  However,  there  are  also  other  Inadequate  intervals  which 
cannot  be  expressed  by  Ad  -  nTj^.  In  experiment  3,  Ad  -  W/x  -  19^. ^36  is  an 
inadequate  interval.  The  relationships  between  Ad  and  Tj^^  are  as  follows; 

Ad  -  0.04  T, 

-  0.16  Tg 

-  2.74  T3 

-  5.48  Tj^ 

-  24.67  T5 

-  49.34  T5 

-  194.44  Tj 

Considering  multiples  of  Ad,  however,  it  becomes  evident 

3  Ad  -  74  T5  -  148  T5 

This  expression  Indicates  that  Ad  and  T^  have  no  common  factor  except  one,  and 
the  observations  d^j,  ^  (d^^^  +  d^jj.^.))  are  uniformly  placed  only  on  three 
different  positions  over  periods  T^  or  T^,  no  matter  how  many  observations  are 
made.  Table  10  shows  the  distribution  of  the  30  observations  used  for 
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obs  4 


Table  10. 


“5 

203.89^ 


“6 

203.894 


obs  # 


“5 

46.691 


-  194.438 

Ad  -  194.. 

0.0021 

0.0044 

0.3387 

0.6774 

0.6785 

0.3498 

0.0126 

0.0252 

0.3494 

0.6990 

0.6841 

0.3683 

0.0157 

0.3554 

0.7110 

0.6843 

0.3688 

0.0216 

0.0432 

0.3551 

0.7104 

0.6906 

0.3813 

0.0185 

0.0371 

0.3529 

0.7061 

0.6853 

0.3709 

0.0177 

0.0356 

0.3538 

0.7078 

0.6837 

0.3677 

0.0181 

0.0365 

0.3515 

0.7033 

0.6837 

0.3677 

0.0183 

0.0369 

0.3570 

0.7144 

0.6916 

0.3836 

0.0243 

0.0489 

0.3619 

0.7241 

0.6975 

0.3953 

0.0325 

0.0655 

0.3726 

0.7455 

0.7063 

0.4131 

“6 

46.691 


0.0662 


0.1254 

0.2508 

0.4247 

0.8494 

0.7219 

0.4438 

0.0169 

0.0338 

0.3154 

O.63O8 

0.61 49 

0.2298 

0.9099 

0.8198 

0.2068 

0.4136 

0.5053 

0.0106 

0.8045 

0.6090 

0.0974 

0.1948 

0.3988 

0.7976 

0.6842 

0.0169 


0.2068 


0.4247 

0.5053 


0.6149 


0.7219 

0.8045 


0.9099 


0.0106 


0.0662 


0.2298 


0.6090 


0.6558 

0.7976 

0.8198 

.8 


An  example  of  the  phase  distribution  over  a  period  produced  by  one 
adequate  interval  of  observations  and  one  inadequate  interval  of 
observations 

-  decimal  part  of  2dj^/Tj^ 
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experiment  3  over  periods  and  Tg  (columns  2  and  3).  It  obvious  that  all 
30  observations  are  approximately  equivalent  to  only  3  observations  to  be  used 
for  solving  the  four  parameters  5A^,  6<>5,  and  6<|(g. 

In  experiment  15  observations  with  intervals  Ad2  -  W2/X  -  ilM.525 
replace  the  sixteenth  to  thirtieth  observations  used  in  experiment  3.  These 
15  new  measurements  were  made  with  W2  -  1*6.691  mgal.  The  relationship  between 
Ad2  and  T^  is  15  Ad2  -  8*1.73  T^.  Table  10  shows  that  the  15  observations  are 
uniformly  distributed  at  15  different  points  over  periods  T5  and  Tg  respec¬ 
tively,  giving  in  total  18  locations  of  observations  over  T5  and  Tg.  The 
accuracy  of  the  results  by  using  both  of  Ad^  and  Ad2  oay  be  seen  to  be  much 
better  than  by  using  only  Adi . 

In  Figure  2  the  phase  distributions  of  observations  over  period  T^  cor¬ 
responding  to  intervals  Adi  -  19*1.1136  and  Ad2  *111.525  are  shown.  From  these 
figures  it  is  clear  that  the  superior  solution  for  the  wave's  parameters  would 
be  the  one  using  measurements  distributed  in  conformance  with  Ad2. 

3  TR(Z  ) 

4.3.2  Solutions  for  the  equation  — r — -  0 

3  Ad 

The  numerical  examples  of  Section  4,3.1  served  to  provide  some  prelirai- 

nau’y  insight  into  the  relationship  between  the  distribution  of  observations 

and  the  accuracies  of  the  adjusted  parameters.  It  is  our  goad,  to  determine 

which  is  the  optimal  uniform  distribution  of  measurements.  To  comply  with  our 

specified  optimization  criterion  we  attempt  to  determine  such  a  Ad  whereby 

TRd^)  -  minimum.  These  Ad  values  have  been  computed  numerically  using  (13) 

and  the  partial  derivative  expressions  (20).  Considering  the  numerical 

procedures  used,  the  computed  roots  are  believed  to  be  in  error  by  less  than 

3  TR(Z  ) 

0.001  c.u.  Several  sets  of  solutions  to  the  equation  _ _  •  0 

3Ad 


40 


are  presented  in  Tables  11a,b.  Notice  each  table  pertains  to  a  different 
periodic  screw  error  function  comprised  of  varying  combinations  of  periods. 

From  Table  11a  it  may  be  seen  that  TR(Zjj)  attains  a  minimum  for  Ad  - 
29.903  c.u.  This  should  be  considered  the  optimum  value  for  the  particular 
periodic  screw  error  function  identified  below  the  table.  Table  11a  also 
illustrates  another  factor  to  consider.  Finding  extreme  points  of  a  function 
using  the  calculus  may  give  a  Ad^  for  which  TR(Ejj)  is  maximum.  As  shown  in 
Table  11b  TR(Zjj)  has  a  minimum  at  Ad  -  I62.8il2  c.u.  for  the  seven  component 
screw  error  function 

Further  information  about  the  relationship  between  the  interval  Ad  and 
the  resulting  minimum  or  maximum  trace  of  appears  in  Table  12.  Some  values 
for  Intervals  Ad^p^^  and  Ad,  were  obtained  from  the  computation  whose  results 
appeared  in  Table  11b.  Table  12  shows  that  although  45  observations  were 
made,  the  equivalent  observations  corresponding  to  Ad,  are  much  less  than  the 
practical  ones;  in  fact,  often  even  less  than  the  minimum  number  of  observa¬ 
tions  required  for  determination  of  unknowns.  With  optimum  intervals  of 
observations,  there  are  almost  no  repeated  positions  of  observation  over  each 
period  Tj^.  All  observations  are  uniformly  or  nonuniformly  distributed  over 
each  period  Involved.  For  example,  a  group  of  observations  with  Ad^pp  - 
162.842  are  uniformly  placed  over  each  period. 

In  a  more  general  setting,  the  relationship  between  Interval  Ad  and 
period  Tj^  may  be  written  in  the  form; 


ROOTS  FOR 

OPTIMUM  Ad 

TR(I^) 

RMS 

AMPLITUDE 

(ugal) 

RMS 

PHASE 

(deg) 

6.590 

0.2992 

1  .3 

18.09 

9.712 

0.2956 

1.3 

17.98 

10.90^ 

0.2929 

1  .3 

17.90 

19.356 

0.2979 

1.3 

18.06 

22.853 

0.2988 

1  .3 

18.08 

26.720 

0.3097 

1.3 

18.41 

27.78^1 

0.2967 

1  .3 

18.01 

29.903 

0.2907 

1.3 

17.83 

35.268 

0.3131 

1  .2 

18.51 

The  following  Ad  give  maximum  TR(Ej^) 

J4.437 

4.0  X  10^ 

49.0 

659.71 

6.700 

1.8  X  10^ 

33.2 

445.0 

8.500 

8.0  X  10^ 

34.8 

435.0 

15.700 

2.9  X  10** 

129.5 

5583.8 

22.900 

7.2 

6.8 

88.48 

24.040 

8.4 

14.8 

96.06 

30.923 

5.0  X  10^ 

52.9 

740.2 

31.900 

6.7 

8.6 

83.53 

35.500 

2.7  X  10® 

18129.5 

5.4  X  10 

Table  11a.  Some  solutions  for  minimum  and  maximum  TR(Ejj).  The  periodic 
screw  error  function  has  the  three  periods 
T^  ■  7.882  C.U.,  T2  ”  3*9^1  c.u.,  T^  •  1.0  c.u. 

W  -  20.0  ragal ,  r  «  21  observations 
Adopt  ”  29.903.  TR(Ejj)  »  minimum 
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TR(Zj^) 


RMS 

AMPLITUDE 

(ugal) 


RMS 

PHASE 

(deg) 


RELATIONSHIP 
BETWEEN  Ad  AND  T,- 


RELATIONSHIP 
BETWEEN  Ad  AND  T. 


n  and  are  mutual  prime  numbers,  i.e.,  n  and  mj^  have  no  common  factor  other 
than  one.  The  sizes  of  n  and  mj^  will  influence  whether  the  interval  Ad  will 
result  in  a  minimum  or  maxirauir  TRiE^^). 

From  this  discussion  it  follows  that  two  sets  of  observations  could  be 
made  with  the  same  weight  but  different  interval  Ad.  At  a  first  glance  they 
may  appear  as  two  completely  different  sets  of  observations,  with  different 
distributions  over  the  range  of  counter  readings.  However,  the  distribution 
of  their  phases  corresponding  to  both  sets  of  observations  over  the  range  of 
any  period  involved  may  be  identical.  This  would  consequently  result  in  the 
same  standard  deviation  of  adjusted  parameters  for  both  sets  of  observations. 

^.3.3  The  effective  observations  and  the  optimum  effective  observations 

We  saw  in  Section  ^.3.2  that  the  accuracy  of  the  parameters  depends  not 
on  the  distribution  of  observations  over  the  range  of  counter  readings. 
Rather,  it  depends  on  the  distribution  of  the  phases  over  the  range  of  each 
period  involved,  provided  the  other  conditions  for  making  observations,  such 
as  the  calibrating  weights,  are  the  same.  The  number  n  in  (41)  defines  the 
number  of  observations  which  occupy  n  separate  phases  over  period  T^.  The  n 
equivalent  observations  under  the  uniform  distribution  are  "effective''  for 
determination  of  unknowns.  Concerning  this,  the  following  definitions  may  be 
establ ished: 

The  effective  observations 

If  the  relationship  between  the  interval  for  a  uniform  distribution  of 
observations  Ad  and  the  period  Tj  can  be  expressed  by  the  form  (41),  the 
number  n  in  the  expression  is  referred  to  as  the  number  of  effective  observa¬ 
tions  for  the  period  T^  . 


A  perfect  group  of  effective  obaervationa 


If  the  number  of  the  effective  observations  are  equal  to  the  practical 
number  of  observations,  namely  n  -  r,  a  perfect  group  of  effective  observa¬ 
tions  for  the  period  is  made. 

The  optimum  group  of  effective  observations 

The  perfect  group  of  effective  observations  for  all  of  the  periods  T^  (i 
1,  2,  k)  is  referred  to  as  the  optimum  group  of  effective  observations, 

and  corresponding  intervals  are  the  optimum  ones. 

A  group  of  observations  with  uniform  distribution  usually  can  not  keep 
their  phase  distribution  uniform  over  all  periods  or  even  one  period,  if  n  •  r. 
There  certainly  are  some  phases  corresponding  to  the  observations  which  nearly 
or  exactly  coincide  with  each  other  over  a  certain  period.  This  means  that 
the  phase  distribution  of  observations  over  the  period  Tj^  would  be  nonuniform 
if  Ad  and  Tj^  can  not  be  nearly  or  exactly  expressed  by  the  form  (41);  i.e.,  a 
pair  of  mutual  prime  numbers  n  and  m^  which  satisfy  (41)  can  not  be  found. 

In  Chapter  3  we  learned  the  criterion  for  optimum  distribution  of  obser¬ 
vations  for  determination  of  the  circular  error  of  one  wavelength  la  to  keep 
observations  uniform  over  the  whole  period.  Let  us  extend  the  criterion  for 
one  wavelength  to  that  for  more  than  one  wavelength.  The  optimum  distribution 
of  observations  for  determining  the  circular  error  with  a  combination  of 
several  wavelengths  la  also  a  uniform  distribution  by  which  all  phases  of 
observations  are  uniformly  distributed  simultaneously  over  each  period  Tj^; 
that  is,  the  interval  of  observations  Ad  should  satisfy: 


The  number  of  observations  r  and  any  number  in  (‘^2)  constitute  a  pair  of 
mutual  prime  numbers. 

Similarly,  the  conditions  for  an  optimum  distribution  over  more  than  one 
wavelength  can  be  written  as  follows 


It  would  be  impossible  to  find  an  interval  Ad  meeting  (43)  if  there  were  no 
common  factors  between  each  period  T^,  for  instance,  between  the  long  wave¬ 
length  -  5325  and  T2  -  1206.  In  practice,  a  better  interval  la  usually 
found  which  is  approximately  fitted  with  (42)  holding  for  most  of  the  periods. 


4.3.4  Design  of  the  optimum  interval  for  uniformly  distributed  observations 
Based  on  the  formula  (42),  the  optimum  interval  of  observations  can  be 
designed.  The  period  T^  is  assumed  to  be  a  long  wavelength  which  is  nearly 
equal  to  the  whole  range  of  counter  readings,  so  it  is  impossible  that  over¬ 


lapping  phases  for  all  observations  occur  over  the  period  Tj^.  Moreover,  the 
period  Tj  is  assumed  to  be  one  mgal  which  is  bound  to  meet  (42).  Hence,  we  do 
not  have  to  take  the  periods  T^  and  Tj  into  account  in  our  design.  It  is  also 
supposed  that  there  exist  some  multiple  relations  between  the  rest  of  the 
periods  Tj^  (i  =■  2,  3.  •••  6); 


T  -  2  T 

‘3  ‘4 
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Using  the  assumptions,  the  optimum  Interval  can  be  determined  by  the  following 


steps: 

1.  Specify  the  number  of  observations,  r. 

2.  Specify  the  range  of  counter  readings,  R,  to  be  calibrated. 


3.  Find  m2^o  -  R/T2. 

4.  Find  an  integer  m2  which  is  nearest  to  m2  g  and  has  no  common  factor 
with  r  except  one. 


5-  Adopt  -  m2  T2/r. 

6.  Specify  Initial  observation  and  then  compile  the  table  of 


(J 


1.  2, 


observations  dof.^.  by  -  ^off.  ,  *  Adopt 

J  j  ' 

This  procedure  may  be  Illustrated  through  the  following  example. 


...  r ) 

Assume  the 


data: 

r  -  44  observations,  T2  -  1206  c.u.,  R  -  6000  c.u.,  02^0  "  >"2  " 

Adopt  -  137.045  C.U.,  doff.  -  100. 3  c.u.,  doff  -  ^off.  ,  *  1  37.045  c.u. 

1  J  J 

(J  -  2,  3.  ...  44),  W  -  180.0  mgal. 

Using  the  above  data,  the  corresponding  standard  deviation  of  the  adjusted 
parameters  were  obtained  by  simulative  computation:  m^^^  -  0.25  x  1 0“^ 
mgal/c .u. 


^2 

"3 

T4 

^5 

^6 

RMS 

”6  A 

0.0033 

0.00070 

0.00031 

0.00062» 

0.00047* 

0.00030 

0.00036 

0.001 3 

0.56 

4.14 

'1. 13 

8.24* 

6.14* 

4.02 

4.64 

5.04 

”y 

0.0032 

0.00069 

0.00031 

0.00061 

0.00047 

0.00030 

0.00035 

0.0012 

0.0030 

0.00069 

0.00031 

0.00062 

0.00046 

0.00336 

0.00035 

0.0012 

(•these  values  of  the  deviations  are  larger  than  the  others  because 


3in  =r-y  and  sin  ^  are  much  less  than  ein  =r-7  (i  "  3i  6,  7)).  It  would  be 

ij^X  T^X 


inconvenient  for  the  gravity  meter's  operator  to  reset  a  counter  reading  from 
the  observation  d^^f.  to  the  next  observation  d^^^  because  Ad^p^  usually  is 
not  equal  or  near  to  ^d^  -  d^^j  -  d^^f  -  W/x.  Therefore,  we  had  better  try  to 
let  the  value  of  Ad^p^.  be  close  to  Ad^. 


4.4  Nonuniform  Distributions  of  Observations 

4.4.1  The  solution  for  the  equations  9TR(Ejj)/3dQff  -  0 

In  order  to  find  the  optimum  nonunifonn  distribution  of  observations  the 
system  of  non-linear  equations: 


(45) 


3  TR(Zj^) 


O 


0 


3  TR(Z^) 


3  d 


off. 


0 


3  TR(Z^) 


d  a 


off 


is  solved  with  an  iterative  method.  Frco  (19)  and  the  expression  (20)  for  the 
partiail  derivatives  with  respect  to  observation  d^f ^ ,  an  expansion  formula  for 


the  left-hand  sides  of  (45)  can  be  obtained.  Several  sets  of  roots  could  be 


found  over  the  given  range  of  each  observation  d^^.^  ,  so  only  those  roots  that 


make  the  trace  minimum  are  taken  as  the  final  solutions  for  (45).  The 


errors  of  the  solutions  are  believed  to  be  less  than  1  c.u. 


since  the  convergence  rate  of  the  iteration  for  solving  the  system  (45) 
is  quite  slow,  it  would  be  very  expensive  to  solve  a  system  which  includes  all 
wavelengths  (k  -  7)  to  be  calibrated.  We  therefore  only  solve  those  systems 
(45)  which  Involve  the  combination  of  less  than  5  wavelengths.  The  simulative 
computations  should  give  us  some  ideas  about  the  distributions.  The  optimum 
nonuniform  distribution  for  determination  of  parameter  sets  x,  and  x, 

are  respectively  investigated.  We  are  also  Interested  in  the  phase 
distribution  corresponding  to  the  optimum  counter  readings  with  weight  off. 

The  phase  distributions  are  illustrated  by  diagrams  similar  in  nature  to  those 
in  Figure  2. 

4.4.2  The  optimum  nonuniform  distribution  for  the  determination  of  parameters 
x ,  ,  and 

Selecting  x,  Y|;,  and  as  the  parameters  for  adjustment,  the  system  of 
equations  (45)  was  solved  by  the  iteration  method.  Some  typical  solutions  for 
the  equations  are  compiled  in  Tables  13a-d.  The  corresponding  optimum  phase 
distributions  are  illustrated  by  the  Figures  3a-d.  In  analyzing  the  results 
the  following  facts  should  be  pointed  out. 

The  trace  of  the  variance-covariance  matrix  of  the  parameters  x,  Yj^,  and 
Zj^  is  reduced  through  the  optimization,  and  the  accuracies  of  almost  all 
parameters  are  improved  to  a  certain  extent.  There  are  no  occasions  when  the 
accuracies  for  some  of  the  parameters  are  improved  while  the  accuracies  of 
others  are  getting  significantly  worse. 

If  the  initial  uniform  distribution  of  observations  for  the  iteration  is 
approximately  equal  to  the  optimum  uniform  distribution,  then  there  are 
practically  no  differences  in  the  accuracies  of  the  parameters  between  the 
optimum  uniform  distribution  and  the  optimum  nonuniform  distribution;  the 
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INITIAL 

FOR 

ITERATION  (c.u. 


180.0  mgal 


STATISTIC 


15  observations 


FOR 

INITIAL  d, 


70.941  c.u 


FOR 

ROOTS 


35.471  c.u 


1000.500 

1007.206 

TR(Exy2) 

2.957  X  10  ^ 

2.847  X  10 
.  -4 

1078.151 

”6x 

4.3  X  10  ^ 

4.3  X  10  ^ 

1133.833 

1208.435 

'"6yl 

0.0005 

0.0005 

1267.167 

1269.299 

”6/2 

0.001 1 

0.001 1 

1 400.500 

1 401 .240 

"'o6y 

0.0008 

0.0008 

1 472.180 

‘°5z1 

0.0006 

0.0005 

1533.833 

1604.709 

'"5z2 

0.001 1 

0.001 1 

1667.167 

1737. 318 

®o6z 

0.0009 

0.0008 

1800.500 

1869.252 

"5AI 

0.0006 

0.0005 

1933.833 

1934.493 

'"5A2 

0.001 1 

0.001 1 

2005.434 

'"06A 

0.0009 

0.0008 

2067.167 

2068.059 

”6<|i1 

8.10 

8.34 

2200.500 

2202.250 

'"64.2 

14.01 

14.14 

2333.833 

2467.167 

2600.500 

2733.833 

2867.167 

2273.189 

2335.061 

2406.002 

2468.046 

2538.046 

2671 .419 

2803.587 

2936.245 

'"06  4. 

11.44 

11.61 

Table  13a.  Some  solutions  for  minimum  TR(Zj^)  taking  x,  as  the  adjusted 

parameters 


1.  • 


INITIAL  OBS 

FOR  ITERATION  ROOTS  FOR  FOR 

(dgff.)  STATISTIC  INITIAL  OBS  ROOTS 


W  =  20.0 

1^-15  T^ 

=■  7.882 

T2  =  3.941 

2000.650 

2002.1 12 

TRdxy^) 

0.1575 

X  10“^  0.2849 

X  10 

2002.650 

2004.106 

3.85  X 

10"^  3.83  X 

10“^ 

2004.650 

2006.174 

®5yi 

0.0006 

0.0005 

2006.650 

2007.842 

'"6y2 

0.0031 

0.001 1 

2008.650 

2009.691 

'"o6y 

0.0022 

0.0008 

2010.650 

201 1 .762 

®6z1 

0.0005 

0.0005 

2012.650 

201 3.893 

”6z2 

0.0023 

0.001 1 

201 4.650 

2015.442 

”o6z 

0.0017 

0.0008 

2016.650 

2016.650 

'"6A1 

0.0006 

0.0005 

2013.650 

2018.705 

”6A2 

0.0013 

0.0J1 1 

2020.650 

2020.760 

"*06  A 

0.0010 

0.0008 

2022.650 

2022.676 

3 

o* 

12.97 

13.19 

2024.650 

2024.650 

"*602 

58.99 

17.17 

2026.650 

2028.650 

2026.650 

2028.672 

"*060 

42.71 

15.31 

INITIAL  OBS 

FOR  ITERATION 

ROOTS 

STATISTIC 

FOR 

INITIAL 

OBS 

FOR 

ROOTS 

W  =  180. 

.0  r  =  21 

T^  =  1206.00 

T,  =  70, 

.941 

T3  =  35.471 

1000.700 

1048.147 

TR ( E^yz  ^ 

0.3566 

X  10“^ 

0.4075  X  10“^ 

1143.557 

1 267.081 

”6x 

40.72  X 

10“^ 

0.36  X  10“^ 

1286.414 

1 390.340 

"’5yi 

0.001 1 

0.0010 

1429.271 

1560.590 

'"5y2 

0.0157 

0.0005 

1572.1 29 

1684.244 

®5y3 

0.0177 

0.0009 

171 4.986 

1812.344 

®o6y 

0.01 37 

0.0008 

1857.843 

191 1 .381 

®6zl 

0.001 1 

0.0010 

2000.700 

2068.074 

“5z2 

0.0486 

0,0005 

2143.557 

2175.790 

“'5z3 

0.0253 

0.0008 

2286.41 4 

2350.749 

®o6z 

0.0316 

0.0008 

2429.271 

2440.683 

®6A1 

0.001 1 

0.0010 

2572.129 

2706.558 

'"5A2 

0.0496 

0.0005 

271 4.986 

2748.960 

®5A3 

0.0296 

0.0009 

2857.843 

2943.760 

®o6A 

0.0334 

0.0008 

3000.700 

3067.046 

®6<1)1 

6.77 

6.08 

3143.557 

3186.691 

”6(1.2 

188.0 

7.06 

3286.41 4 

3287.745 

“6(1.3 

118.0 

11.97 

3429.271 

3443.746 

”o6(1. 

128.2 

8.76 

3572.129 

3690.705 

371 4.986 

371 4.986 

3857.843 

3994.513 

Table  13c.  Some  solutions  for  minimum  TR(Ejj)  taking  x,  as  the  adjusted 

parameters 


INITIAL  OBS 
FOR  ITERATION 


ROOTS 


STATISTIC 


FOR 

INITIAL  OBS 


FOR 

ROOTS 


W  -  180.0 


T,  =  1206.00 


70.9^1 


35.471 


100.300 

226.884 

TR(E^y^) 

0.1804  X  10"^ 

0.4073  X  10"^ 

386.01 4 

525.737 

mfiy 

1.773  X  lO"^ 

0.360  X  10'^ 

671 .729 

785.606 

"“Syi 

0.0010 

0.0010 

957.443 

1020.733 

®6y2 

0.0016 

0.0005 

1243.157 

1311.965 

“6y3 

0.0026 

0.0009 

1528.871 

1596.750 

“oay 

0.0019 

0.0008 

181 4.586 

2075.035 

™6z1 

0.0010 

0.0010 

2100.300 

2308.249 

"'6z2 

0.0021 

0.0005 

2386.01 4 

2490.227 

"’5Z3 

0.0015 

0.0009 

2671 .729 

2915.239 

®o6z 

0.0016 

0.0008 

2957.443 

3024.489 

"eAi 

0.0010 

0.0010 

3243.157 

3310.101 

'"5A2 

0.0025 

0.0005 

3528.871 

3667.872 

""(SAS 

0.0024 

0.0009 

381 4.586 

3850.079 

'"o6A 

0.0021 

0.0008 

4100.300 

4166.964 

'”64'1 

6.24 

6.06 

4386.014 

4420.809 

^602 

10.53 

7.05 

4617.729 

4739.932 

^603 

22.93 

1 1  .91 

4957.443 

5026.791 

"^060 

15.01 

8.72 

5243.157 

5243.157 

5528.871 

5598.793 

5814.586 

5888.1 34 

Table  13d. 

Some  solutions  for 
parameters 

minimum  TR(j: 

taking  x,  , 

as  the  adjusted 
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change  in  phase  distribution  is  also  not  significant  (Tables  13ai  13c,  and 
Figures  3a.  3c). 

If  the  initial  uniform  distribution  of  observations  for  the  iteration  is 
far  from  the  optimum  uniform  distribution,  namely,  the  corresponding  interval 
of  observations  is  inadequate,  then  the  accuracies  of  all  parameters  will  be 
considerably  improved  by  optimization  with  the  optimum  nonuniforra  distribution 
of  observations.  The  dense  phase  distribution  located  at  a  few  positions  over 
the  periods  corresponding  to  the  Inadequate  intervals  will  nearly  uniformly 
scatter  over  all  periods  and  change  into  an  optimum  nonuniform  distribution 
which  is  similar  to  the  uniform  distribution  (of  Tables  13b,  1  3cl  and  Figures 
3b,  3d).  According  to  the  above  facts,  we  would  consider  that  the  accuracies 
of  the  parameters  x,  ,  and  obtained  from  optimization  with  the  uniform 
distribution  of  observations  can  hardly  be  improved  by  optimization  with  the 
nonuniform  distribution. 

It.  4. 3  The  character  of  the  optimum  nonuniforra  distribution  for  determining 
the  parameters  x,  6A^  and 

Tables  1 4a-c  and  Figures  4a-c  show  the  comparison  in  phase  distribution 
and  accuracy  between  two  types  of  optimum  distributions  of  observations. 
Several  facts  are  evident.  First  the  trace  of  the  variance-covariance  matrix 
of  the  parameters  x,  5Aj^  and  forraed  from  the  adjustment  of  nonuniformly 
distributed  observations  in  each  case  is  markedly  less  than  with  the  optimum 
uniform  distribution.  Secondly,  the  accuracy  of  the  unknown  phase  lag 
correction  is  Improved  by  the  adjustment  with  the  optimum  nonuniforra 
distribution  of  observations.  The  standard  deviation  m^^  is  reduced  about  0.5 
-  3.0  degrees  for  these  examples.  On  the  other  hand  a  significant  reduction 
f  the  accuracies  of  the  other  parameters  5Xj^,  6Aj^,  and  occurred.  The 
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INITIAL  OBS 

FOR  ITERATION 

ROOTS 

FOR 

FOR 

^“^off  ^ 

^'^off  ^ 

STATISTIC 

INITIAL  OBS 

ROOTS 

(c.u. ) 

(c.u. ) 

W  =  180.0  mgal  r  -  15  observations  =  5325.00  c.u.  T2  »  1206  c.u. 


300.300 

129.914 

TR(xA4,) 

0.01627 

0.01072 

700.300 

726.907 

“ax 

0,43  X  10"^ 

0.52  X  10"^ 

1100.300 

»  900.300 

'''5A1 

0.0054 

0.0054 

1500.300 

1333.133 

™6A2 

0.0012 

0.0026 

1900.300 

1936.407 

^06  k 

0.0039 

0.0043 

2300.300 

*2100.300 

“a^i 

0.96 

0.95 

2700.300 

2535.219 

®<S<t)2 

7.24 

5.76 

3100.300 

3140,650 

"'o6(|» 

5.16 

'J.13 

3500.300 

*3300.300 

®y1 

0.0052 

0.0053 

3900.300 

3740.352 

®y2 

0.0012 

0.0012 

4300.300 

4344.585 

""oy 

0.0038 

0.0038 

4700.300 

4500.300 

™z1 

0.0051 

0.0050 

5100.300 

4951 .614 

™22 

0.0012 

0.0025 

5500,300 

5549.644 

"oz 

0.0037 

0.0040 

5900.300 

*5700.300 

Table  14a. 

Some  solutions  for 

minimum  trace 

(E^)  taking  x,  6A 

^  and  as 

adjusted  parameters. 

3TR(i:  ) 

(*  not  a  root  of 

3d  °  ' 

this  observation 

is  optimum  : 

^  off 

the  given  range  for  finding  optimum  points.) 


INITIAL  OBS 


FOR  ITERATION 
^^off ^ 

ROOTS 
(^off ^ 

STATISTIC 

FOR 

INITIAL  OBS 

W  =  180.0 

r  »  15  T,  = 

70.941  T2  =  35, 

.472 

2034.033 

2027.071 

0.076840 

0.05153 

2100.700 

21 18.630 

'"6X 

0.44  X  10"^ 

0.64  X  10“^ 

2167.367 

2134.658 

''6A1 

0.00058 

0.00062 

2234.033 

2225.730 

'"6A2 

0.00116 

0.00351 

2300.700 

*2267.367 

”o6A 

0.00092 

0.00252 

2367.367 

2382.799 

®6(J)1 

8.03 

7.53 

2434.033 

2438.405 

®6(J)2 

13.70 

10.60 

2500.700 

2509.002 

™06()) 

1  1  .23 

9.19 

2567.367 

2580.209 

®y1 

0.00051 

0.00051 

2634.033 

2631.732 

®y2 

0.00106 

0.00254 

2700.700 

2702.695 

0 

a 

0.00083 

0.00183 

2767.367 

2736.943 

®Z1 

0.00058 

0.00060 

2834.033 

2827.802 

®z2 

0.00114 

0.00256 

2900.700 

2967.367 

2878.777 

2969.362 

®oz 

0.00090 

0.00186 

Table  l^b.  Some  solutions  for  minimum  trace  (Ej^)  taking  x,  6A^  and  as 
adjusted  parameters. 

3TR(E^) 

(*not  a  root  of  — r-r -  •=  0  ;  this  observation  is  optimum  in 

^  off 


the  given  range  for  finding  optimum  roots.) 


INITIAL  OBS 

FOR  ITERATION 

ROOTS 

FOR 

FOR 

(doff) 

^^off ^ 

STATISTIC 

INITIAL  OBS 

ROOTS 

W  =  180.0 

r  =  15  T-j 

-  7.882  T2  =  3 

.941 

2001 .650 

*2000.650 

0.2070 

0.0960 

2003.650 

2003.785 

'"5X 

3.85  X  10"^ 

5.44  X  10”^ 

2005.650 

2005.326 

®6A1 

0.00053 

0.00063 

2007.650 

2007.509 

®6A2 

0.00364 

0.00278 

2009.650 

2009.310 

®o6A 

0.00260 

0.00201 

201 1 .650 

201 1 .574 

‘"6i>1 

13.54 

1 1.75 

2013.650 

2013.115 

%<t>2 

22.28 

13.30 

2015.650 

2015.481 

18.44 

12.55 

2017.650 

2017.249 

myi 

0.00053 

0.00063 

2019.650 

2019.578 

®y2 

0.00239 

0.00114 

2021  .650 

2021 .188 

®oy 

0.00173 

0.00092 

2023.650 

2023.529 

®Zl 

0.00055 

0.00048 

2025.650 

2025.326 

3 

N 

0.00307 

0.00286 

2027.650 

2027.578 

®02 

0.00221 

0.00191 

2029.650 

2029.318 

Table  14c. 


Some  solutions  for  minimum  trace  taking  x,  6A^  and  as  the 

adjusted  parameters. 

3TR(E^) 

(*not  a  root  of  — rr - —  =  0  ;  this  observation  is  optimum  in 

^'^off 

the  given  range  for  finding  optimum  roots.) 


.  UMIFARM 
aiSTRItUTION 


Fi Kure 


4a.  Ti  =  5325.00  c.u.  r  =  15  observations  Refer  Table  14a 


4a.  T2  =  1206.00  c.u.  r  =  15  Refer  Table  14a 


4.  The  contrasts  between  the  phase  distribution  corresponding  to 
unifoni  distribution  of  observations  and  optimal  nonunifom 
distribution  of  observations,  taking  x,  6A,  64  as  paraneters. 


4b.  T,  =  70.941  c.u.  r  =  15  Refer  Table  14b 


35.471 


Hefor  Table  14b 


losses  in  the  accuracies  of  the  parameters  SA^,  and  reach  about  0.5  - 

2.5  ugal ,  and  the  linear  scale  factor  5x  about  0.2  x  10~^  -  1.6  x  10~^ 
c.u./mgal.  This  means  that  it  is  advantageous  only  for  the  determination  of 
the  phase  lag  to  solve  the  system  (45).  However,  it  is  Just  in  the  sense 
of  comparing  optimally  distributed  uniform  and  nonuniform  observations. 

In  order  to  realize  this  situation  mentioned  above,  let  us  analyze  the 
character  of  the  optimum  nonuniform  distribution  of  observations.  In  Figs. 
4a-c,  there  is  a  common  character  to  the  optimum  nonuniform  distribution. 

That  is,  most  of  the  observations  are  concentrated  near  -  0“  and  180°  over 
each  period  with  the  exception  of  the  long  wavelengths,  which  resembled  the 
phase  distribution  of  the  uniform  distribution. 

The  phase  distribution  for  wavelengths  T^  and  T^  are  admost  the  same. 

The  mean  positions  of  some  dense  distributions  represented  by  average  phase 
angle  9„,  are  approximately  equal  to  45°,  150°,  225“  and  330°.  For 
wavelengths  T^  and  T^,  -  15*  and  165“.  Regarding  wavelength  T2,  the 

corresponding  phase  distribution  also  tends  toward  0°  and  180°,  but  not 
strongly. 

In  a  typical  observation  equation,  the  coefficient  for  has  a  factor 

which  is  usually  much  less  than  1  mgal.  It  is  obvious  that  the  accuracy  of 

determination  of  would  be  Improved  if  the  value  of  cos(^  d^  -  (jj.)  were 

1 

increased:  l.e.  if  Q  -  d  -  (ji,  •  0°  or  180°.  Otherwise,  the  coefficient 

m  T  m  1 

for  6(>j^  would  be  small  and  the  accuracy  of  the  adjusted  would  correspond¬ 
ingly  decrease.  This  implies  that  the  variance  of  is  more  sensitive  to  the 
phase  Gjjj  than  the  variance  5A.  Concerning  wavelength  T^  ,  the  corresponding 
phase  distribution  is  nearly  uniform.  There  is  no  significant  difference  in 
the  accuracy  of  5A  and  5i^,  and  especially  6<>,  between  the  optimum  nonuniform 
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and  uniform  distributions.  As  discussed  previously,  a  uniform  distribution  of 


observations  is  preferable  for  determination  of  scale  factor  x  and  amplitude 

Ai. 

It  was  pointed  out  in  section  3.2  that  the  condition  sinZGj^j  -  0  ( j  »  1  , 
2,  ...  r)  must  be  Imposed  for  the  simultaneous  minimization  of  TR(6Aj^)  and 
TR(6(>j^).  However,  this  type  of  conditional  distribution  could  not  be 
satisfied  for  all  wavelengths  in  the  periodic  error  function  for  either  a 
nonuniform  or  uniform  distribution.  This  would  be  one  of  the  reasons  why  the 
simultaneous  minimization  of  TR{5Aj^)  and  TR(5(J)£)  could  not  generally  be 
reached  in  the  optimization  procedure,  especially  in  the  ones  with  a 
nonuniform  distribution  of  observations.  For  the  uniform  distribution, 
however,  a  simultaneous  minimization  of  TR(SAj^)  and  TR(6(>j^)  is  nearly 
realized,  since  both  traces  decrease  at  comparable  rates  Ih  the  optimization 
procedure. 

In  addition,  it  is  very  important  that  the  maximum  effect  of  m^^  on 
converting  counter  readings  to  their  milligal  equivalent  Is  itself 
However,  the  error  of  (p,  causes  the  maxlmun  error  of  the  milligal 

equivalent  to  be  A  for  example,  if  -  3°  -  0.05  (rad)  and  A  -  10 

ugal ,  then  A  m^^  -  0.5  ugal.  As  Is  mentioned  previously,  a  3  degree  decrease 
of  m^^  would  result  in  a  3  ugal  Increase  of  m^^  in  the  optimization  procedure 
with  the  nonuniforra  distribution.  In  view  of  this  point,  it  seems  impossible, 
in  practice,  to  obtain  an  optimum  distribution  for  highly  accurate  adjusted 
parameters,  if  only  the  optimization  procedure  under  a  nonuniform  distribution 
of  observations  is  employed. 


4.5  Requisite  Number  of  Observations 

Usually,  a  complete  procedure  of  calibration  expends  much  time,  for 
instance,  2-3  weeks,  and  therefore  involves  appreciable  expense.  We  should 
consider  the  important  question,  "how  many  observations  are  needed  for  deter¬ 
mining  the  circular  error  with  a  reasonable  accuracy?"  With  respect  to  one 
wavelength,  an  answer  has  been  given  to  this  question  in  Section  3.3. 
Concerning  the  combination  of  more  than  one  wavelength,  it  is  difficult  to 
answer  the  question  exactly  and  absolutely.  However,  an  approximate  estimate 
of  the  requisite  number  of  observations  can  be  formed  by  some  simulative 
computations . 

First  the  accuracy  criterion  has  to  be  given  for  the  estimation  of  the 
number  of  observations  used  in  the  simulative  computation.  The  accuracy  of 
determination  of  amplitude  for  the  calibration  of  an  LCR  "G"  gravity  meter  has 
been  in  the  range  of  about  0.4  ygal  to  6.5  ygal ,  and  2.0  usal  in  mean  value. 
The  accuracy  of  determination  of  phase  lag  in  the  "G"  gravity  meter  has  ranged 
from  3®  to  29°,  14®  In  the  mean  value  (Kanngieser  and  Torge,  1981).  As  is 
known,  the  minimum  division  of  the  dial  in  a  "G"  gravity  meter  is  roughly 
equivalent  to  10  ugal.  The  maximum  error  of  the  readings  can  be  considered  to 


be  0.005  G.u.  It  would  be  preferable  that  the  RMS  m^^^  should  be  estimated  by 
the  RMS  of  Involved  in  the  periodic  screw  error  function.  According  to  the 
empirical  data  from  the  calibration  for  LCR  "G"  gravity  meters  (Kanngieser  and 


Torge,  1981),  the  range  of  the  amplitudes  have  been  found  which  are  as 


follows : 


<  5  ugal  for  the  1  c.u.  period 

<  2  ugal  for  the  3.9  c.u.  period 

<  9  ugal  for  the  7.9  c.u.  period 

<  3  ugal  for  the  35  c.u.  period 


<  20  ligal  for  the  71  c.u,  period 


•i'. 

"  •  •  • 


<  20* 

(*  found  using 

'  M  s' 

10.7  ugal.  If 

follows : 

<11° 

b 

<  29® 

my. 


A'..', 


m 

m 


N*.  % 


“  i  • 


SV-' 

A' 


A- 

JL.'  t. 


<•- 

>y-. 

V'. 


MV 


p: 

tKv 


m 

mi 


(i-v-i 


are  as 


<  6°  for  the  7.9  c.u.  period 

<  7°  for  the  35  c.u.  period 

<  3°  for  the  71  c.u.  period 

<  3“  for  the  1206  c.u.  period 

The  root  mean  square  of  m^^^  is  required  to  be  less  than  13®.  From  the  ampli¬ 
tudes  given  for  simulative  computation  in  this  report,  the  limiting  m^^  is 


15®,  i.e,,  '  15®. 


Figures  5a-h  show  some  relationships  between  the  number  of  observations 
and  the  RMS  values  of  the  amplitudes  and  phase  lags,  as  well  as  the  maximum 
RMS  in  the  RMS  of  the  adjusted  quantity  Ad  -  d^^^  -  d^ff.  The  curves  in  Figure 
5a  are  drawn  with  the  data  from  the  simulative  computation  under  current 
calibration  procedures  and  apparatus.  30  observations  using  the  180  mgal 
weight  are  made  at  first,  cind  then  we  increase  sequentially  the  observations 
using  the  20  mgal  weight  at  intervals  of  2  observations  up  to  more  than  300. 

It  is  found  that  a  decrease  of  m^^  from  1.il3  to  1.32  ugal  occurs  along  with 
the  increase  of  the  number  of  observations  from  up  to  32^.  In  other  words, 
an  Increase  of  281*  observations  brings  about  an  improvement  in  the  accuracy  of 
adjusted  parameters  only  by  0.11  ugal,  and  m^^^  has  remained  over  1.3 
ugal.  The  Increase  of  observations  using  only  a  small  weight,  such  as  20 
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Figure  5.  The  relationship  between  the  standard  deviation  of  adjusted 
paraaeters  and  the  nuaber  of  observations. 


mgal ,  could  not  lead  to  a  significant  improvement  in  the  accuracy  of  parameter 
6Aj^.  However,  the  improvement  in  the  accuracy  appears  remarkable  as  the 
observations  using  a  large  weight,  for  instance,  180  mgal,  are  increased. 

Let  us  have  a  look  at  Figs.  5a-h,  the  figures  indicate; 
r  >  60,  m^^  <  1.3  pgal 

for  a  weight  of  180.0  mgal 

r  >  90,  m^^  <  1.0  ygal 

r  >  40,  m^^  <  1 .3  pgal 

for  a  weight  of  178.962  mgal 

r  >  66,  m^^  <  1.0  ugal 

r  >  30,  m^^  <  1.2  ugal 

for  a  weight  of  242.672  mgal 

r  >  45 ,  m^^  <  1.0  ugal 

If  <  1.3  pgal  is  required,  about  40  to  60  observations,  for  which  a  large 
weight  should  be  adopted,  are  needed.  If  m^j^  <  1.0  ugal,  about  70  to  90 
observations  are  needed.  In  this  case,  the  maximum  RMS  of  the  adjusted 
milligal  equivalent  to  the  counter  readings  difference  dQ^j  -  would  be 

less  than  about  1.8  ugal.  If  a  larger  weight,  for  instance,  242.672  mgal  were 
used  Instead  of  the  180.0  mgal  weight,  the  number  of  required  observations 
could  be  less  than  45  to  meet  the  mentioned  accuracy.  From  Figures  5  we  note 
that  about  80  to  90  observations  would  be  enough  when  m^^  <  1.3  ugal  is 
specified.  This  requires  30  extra  observations  using  a  weight  larger  than  180 
mgal,  and  50  to  60  extra  ones  using  a  smaller  weight  of  about  50  to  80  mgal. 
According  to  the  simulative  computations  in  Sections  4.2  and  4.3,  the 
requisite  number  of  observations  would  be  less  than  the  above  estimated  number 
of  observations  required,  provided  the  weight  and  the  distribution  of 
observations  are  selected  and  designed  by  the  criteria  of  optimization. 
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Concerning  the  requleite  observations  for  determining  adjusted  parameter 
phase  lag  with  the  reasonable  accuracy  of  say  m,^  <  130,  it  is  obvious,  from 
the  figures  that  the  accuracy  would  be  easily  obtained  provided  about  30  to  AO 
observations  are  made  with  a  proper  Interval  for  uniform  distribution  and 
using  an  appropriately  large  weight  in  the  calibration.  We  emphasize  that  the 
intervals  and  weights  have  to  be  carefully  selected  based  on  those  criteria 
that  were  discussed  in  Sections  4.2  and  4.3.  In  addition,  the  -  r  curves 
in  the  figures  show  that  the  values  of  mj^j,  are  basically  inversely 
proportional  to  the  number  of  observations. 

Finally,  we  point  out  that  the  results  of  this  section  just  come  from 
simulative  computations  with  a  specific  model  for  the  periodic  screw  error 
function,  where  it  is  assumed  that  the  standard  deviation  of  an  observation  is 
equal  to  0.002  c.u.  If  there  would  be  some  change  in  these  conditions  for 
practical  calibration,  the  results  might  have  some  small  differences. 

5.  CONCLUSION 

The  sizes  of  the  weights  used  in  calibration  have  significant  effects  on 
the  accuracies  of  all  adjusted  parameters  of  the  periodic  screw  error 
function,  since  the  standard  errors  of  the  parameters  are  sensitive  to  the 
quantity  sln(TTW/T^x ) .  The  optimum  weights  for  detecting  the  periodic  errors 
should  have  the  sizes  of  nearly  (n^  +  -i]  Tj^x;  namely,  the  value  of 
sin(iT  W/T^x)  should  be  approximately  one.  The  larger  the  value  is,  the  better 
the  corresponding  weight.  In  the  case  of  calibrating  the  periodic  screw  error 
function  with  one  wavelength,  the  optimum  size  of  the  weight  is  equal  to 
(n  +  •^)  Tx,  no  matter  what  type  of  parameters  are  selected  in  the 


adjustment.  For  several  wavelengths,  it  is  impossible  to  find  only  one  weight 

having  the  special  size  for  fitting  each  period  involved.  As  a  general 

rule,  the  optimum  weights  could  be  found  by  solving  the  equation  3TR(E^)/3W  =  0 

for  minimizing  TRCZj^).  The  study  in  this  report  shows  that  the  solutions  for 

the  equation  are  approximately  equivalent  to  the  solutions  for  the  equation 
k  IT  W 

3[  Z  sin  -= — =  0.  For  two  wavelengths  fulfilling  Ti  -  ZToi  the  roots 
i  =  1  i 

of  both  equations  are  expressed  by  -  (n  +  0.32)  x  .  For  more  than  two 

wavelengths,  in  order  to  obtain  higher  and  homogeneous  accuracy,  it  is 
preferable  that  one  weight  and  another  compensative  weight  W2  are  adopted 
so  that  the  value  [  sin  (tt  W^/Tj^x)|  +  |  sin  (u  W^/T^x)!  is  greater  than  1.6 
(see  Tables  10  and  11),  except  for  long  wavelengths  (i  »  1,  2), 

Since  the  selection  for  the  weights  is  related  to  the  linear  scale  factor 
X  and  there  are  some  differences  in  scale  factors  between  different  gravity 
meters,  a  meter  should  have  its  own  weights  as  an  accessory  of  the  meter.  In 
addition,  a  large  weight  can  be  also  used  for  detecting  a  very  short  periodic 
screw  error  provided  the  corresponding  value  |  sin  (itW/Tx)|  is  great 
enough.  A  large  weight  is  better  than  a  small  weight  for  improving  the 
accuracy  of  linear  scale  factor  x  and  amplitude  A.  If  a  non-linear  term  of 
the  scale  factor,  such  as  quadratic  term,  can  be  neglected,  it  would  be  appro¬ 
priate  to  employ  one  or  two  large  weights  for  the  calibration. 

The  optimum  distribution  of  observations  for  calibrating  the  periodic 
screw  error  with  one  wavelength  is  a  uniform  distribution.  If  Y  -  A  coscji  and 
Z  »  A  3in(>  are  taken  as  parameters  in  the  adjustment,  the  optimum  uniform 

distribution  has  to  satisfy  the  condition  expressed  by  the  equations 
r  r 

Z  CO320  «  E  3in20  =■  0  [0  '  r  (d  +  d  )]  over  the  range  of 

,  ,  om ,  .  ,  om ,  om ,  T  of  f .  on . 

J-1  J  J-1  J  J  J  J 

one  period.  In  this  case,  the  highest  accuracy  for  determining  the 


parameter's  amplitude  A  and  phase  lag  <p  would  not  be  attained,  especially  when 
A  is  very  small.  However,  the  difference  in  accuracy  between  the  two  sets  of 
unknown  parameters  is  not  significant. 

For  determining  the  parameters  A  and  <f>,  the  optimum  distribution  of 

observations  is  required  to  satisfy  the  stronger  conditions  3in20j|,j  •  0 

(0^  -  i  ^  ~  minimizing  both  of  TRdg;^)  and  TRCZg^).  The 

condition  means  that  the  distribution  of  observation  is  defined  by  d_  » 

J 

(doff  ^  '^on.)^2  -  (n/i4  >  ^/2-n)  T. 

J  J 

With  regard  to  the  optimum  distribution  of  observations  for  detecting  the 
periodic  errors  with  more  than  one  wavelength,  two  types  of  distributions, 
uniform  and  nonuniform,  are  respectively  discussed.  For  the  uniform  distribu¬ 
tion,  the  Interval  of  observations  should  fit  in  with  all  periods  appearing  in 
the  periodic  error  function,  especially  with  some  middle  and  short  wave¬ 
lengths.  Otherwise  the  number  of  effective  observations  would  be  much  less 
than  the  practical  number  of  observations.  As  a  matter  a  fact,  it  is  the 
phase  distribution  over  the  range  of  periods  involved  and  corresponding  to 
observations  that  acts  in  determining  the  standard  errors  of  the  adjusted 
parameters  and  not  the  distribution  of  observations  over  the  range  of  counter 
readings.  The  uniform  distribution  of  observations  could  result  in  nonuniform 
phase  distribution.  The  number  of  separate  phases,  which  is  referred  to  as 
the  number  of  effective  observations,  is  always  less  than  the  number  of 
practical  observations,  or  at  most  equal  to  them.  The  number  of  the  phases  n 
distributed  over  a  period  Tj^  can  be  found  by  the  expression  n  Ad  -  M.  •  T^  (n 
and  are  two  mutual  prime  numbers).  If  the  nixnber  n  is  much  less  than  the 
number  of  observations  r,  then  the  corresponding  interval  Ad  is  considered  to 
be  inadequate  for  the  distribution.  When  the  nimber  n  is  just  equal  to  the 


number  r  for  all  periods  occurring  in  the  periodic  error  function,  the  corre¬ 
sponding  interval  Ad  is  optimum,  so  is  the  corresponding  distribution.  The 
optimun  Adopt  can  be  determined  if  sane  multiple  relation  exists  between  the 
periods  (see  Section  4.3.J4).  For  LCR  "G"  gravity  meter  having  the  old  gear 
box,  relations  between  the  periods  of  1206,  1206/17,  1206/34,  1206/153, 
1206/306  and  1  are  present.  All  observations  are  always  effective  for  the 
very  long  period  being  equal  nearly  to  the  total  range  of  7000  mgal .  The 
optimum  intervals  for  uniform  distribution  can  also  be  found  by  solving  the 
equation  3TR(Zjj)/3Ad  -  0.  The  solutions  for  the  equations  are  basically 
consistent  with  those  intervals  that  are  found  from  the  expression  r  Ad  «  T 

With  regard  to  the  nonuniform  distribution  of  observations,  the 
simulative  computations  for  optimization  indicate  that  the  optimum  nonuniform 
distribution  of  observations  does  not  have  more  advantages  for  an  improvement 
In  accuracy  of  adjusted  parameters  than  the  optimum  uniform  distribution.  If 
it  is  the  trace  of  the  covariance  matrix  of  the  parameters  x,  Y  and  Z  that  is 
minimized  in  optimization,  there  is  no  significant  difference  in  accuracies  of 
the  parameters  between  the  optimum  nonuniform  distribution  of  observations  and 
the  optimum  uniform  distribution. 

Minimization  of  the  trace  of  covariance  matrix  of  the  parameters  x,  5A 
and  61^  in  optimization  is  not  preferable  for  improving  accuracy  of  all  para¬ 
meters  entering  into  the  adjustment.  Although  the  accuracy  of  the  parameter 
5i^  can  be  improved  significantly  in  the  optimization,  the  accuracies  of  the 
parameters  x  and  6A  would  become  worse,  and  even  below  the  requirement  of  the 
accuracy  for  the  calibration.  The  reason  is  that  the  trace  i-3  more 

sensitive  to  the  phase  distribution  of  observations  than  trace  in  the 

process  of  optimization  by  solving  to  the  equation  3TR  ( ) /3  dg^  -  0.  The 


phase  angles  corresponding  to  the  observations  would  be  concentrated  on  some 
special  phase  positions  over  a  period,  such  as  0°,  180°,  etc.  The  dense  phase 
distribution  results  in  reducing  the  accuracy  of  the  parameters  x  and  6A  which 
basically  represent  the  accuracy  of  the  conversion  from  counter  readings  into 
milligal  equivalents. 

This  leads  to  the  conclusion  that  so  called  optimum  nonuniform 
distribution  of  observations  would  not  be  preferable  for  the  calibration  to 
the  optimum  uniform  ones,  in  the  sense  of  optimization  criterion  discussed 
previously . 

How  many  observations  are  needed  for  the  calibration  of  periodic  screw 
errors  with  reasonable  accuracy?  For  detecting  the  error  with  one  wavelength, 
the  study  indicates  that  the  TR(£^^)  and  TRCE^^)  are  both  inversely  propor¬ 
tional  to  the  number  of  observations  r  made  with  optimum  distribution  and 
using  optimum  calibrating  weights.  If  the  RMS  error  of  5A,  m^^^,  is  required 
to  be  less  than  0.001  mgal ,  and  o  -  0.002  c.u.  is  assumed,  then  the  number  of 
observations  made  in  the  calibration  has  to  be  more  than  5  in  theory.  In 
general,  letting  m^^^  -  A  then  the  requisite  number  of  observations  r  for 

calibrating  the  periodic  error  with  one  wavelength  can  be  approximately  esti- 

2 

mated  by  the  formula  r  >  (xo/(3in  A9c))  ,  where  e  is  a  specified  accuracy  for 

adjusted  parameter  A  and  A0  -  irW/T, 

Concerning  the  problem  about  necessary  observations  for  calibrating  the 
errors  with  more  than  one  wavelength,  the  estimate  is  made  by  the  simulative 
studies.  Based  on  the  error  of  a  counter  reading  of  an  LCR  "G"  gravity  meter 
and  the  empirical  data  of  standard  errors  of  adjusted  parameters  documented  by 
a  great  deal  of  practical  calibrations  for  LCR  "G"  gravity  meter  (Kanngieser 
and  Torge  1981,  and  Becker,  1981),  it  would  be  appropriate  that  the  error  m^^ 


or  i3  required  to  be  leas  than  5//^  ugal .  where  k  la  the  number  of  the 


periods  involved  in  the  periodic  errors.  Assuming  k  -  7,  then  m^^  and 


i 


should  be  less  than  1.3  ugal  or  1.0  pgal ,  and  m^j^^  less  than  13°.  which  are 


considered  as  an  accuracy  criterion  for  the  estimate. 

The  simulative  studies  indicate  that  when  a  large  weight,  such  as  180.0 
mgal ,  is  used  in  the  calibration,  the  RMS  errors  of  all  adjusted  parameters  x, 
6A  and  5<^  rapidly  decrease  with  the  increase  of  observations.  But  the 
increase  of  observations  only  using  a  small  weight,  e.g.  20.0  mgal,  can  hardly 
result  in  an  Improvement  in  the  accuracy  of  parameters  6A  and  x.  For  in¬ 
stance,  an  Increase  of  28U  observations  from  1*0  to  324  including  30  observa¬ 
tions  using  weight  180.0  mgal,  brings  about  an  improvement  in  RMS  error  of  6A 
only  by  0.14  ugal  and  the  RMS  error  of  6A  has  remained  more  than  1.3  ugal,  and 
RMS  error  of  linear  scale  factor  x  keeps  near  a  constant  0.03  x  10~^  mgal/c.u. 
(a  -  0.02  mgal ) . 

If  a  non-linear  term  of  scale  factor  is  negligible,  those  observations 
that  are  made  with  a  large  weight  are  considered  for  estimating  the  requisite 
nvinber  of  observations.  By  the  simulative  studies,  about  40  to  60  observa¬ 
tions  using  nearly  180  mgal  weight  are  needed  for  keeping  the  RMS  error  of  6A 
less  than  1.3  ygal ,  and  about  70  to  90  observations  for  keeping  it  less  than 
1.0  ugal.  If  a  weight  larger  than  l80  mgal,  such  as  242.672  mgal  is  adopted, 
about  35  observations  are  enough  for  having  m^^  less  than  1.3  ugal,  and  about 
45  observations  for  having  less  than  1.0  ugal. 

The  requisite  observations  could  be  about  5  to  10  less  than  the  above 
estimated  number  of  observations  required  when  the  weight  and  the  intervals 
for  uniform  distribution  of  observations  are  strictly  selected  by  the  criteria 


for  optimization  presented  in  the  report. 
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APPENDIX  A 


Derivation  of  the  formulae  for  the  partial  derivatives  of  Adg,  AC  and  AS  with 
respect  to  Ad  and  W 


The  function  for  converting  a  counter  reading  d  to  its  milligal 
equivalent  D  is: 


D  =  xd  +  EAcos  (Xd  -  41)  =  f(d) 


The  inverse  of  the  function  f(d)  is  d  =•  F(D) 

f'(d)  =  ^  =  X  -  ZXAsin  (Xd  -  4,) 
dd 

Let  m  -  f ' (d) ,  (or  ffloff)  =  f  (dQf.f.),  m2  (or  m^^)  =  f'(dQ„) 


^on  =  -  F(f(doff.)  *  W)  =  F(Doff.  *  W) 

‘^'^0  ”  '^on  ~  '^off 

■  '^Offj 

The  above  expressions  are  the  basic  formulae  for  the  derivation  below 


3Ad. 


'^^on  -  Vf' 
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off 
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APPENDIX  B 


Derivation  of  the  formulae  for  the  Trace  (Zj^)  formed  by  the  adjustment  for 
determining  the  parameters  of  a  periodic  screw  error  function  with  one 
wavelength 


The  formulae  for  the  Trace  (.1^)  is: 


TR(Z  )  -  — 
X  r 


2o^  [  I  43^  •  t  4C^  I 

°  J-1 


r  r 

Z  AS^  Z  AC^  -  [  f  AZS  AC  >12 
J  J=1  J  J=^  J 


o  2  2 

2o  X  r 
o 


r  r 

sin^  AG  [r^  -  (  Z  cos20_  [  E  3in20_  )^] 

.  ,  om  '  \  ,  om . 

J“1  J  J-1  J 


expand  the  function  f(d)  (see  appendix  A) 


f(d)  =  D  =  X  d  +  A  cos(y^  d  -  <>] 


«  X  d  +  Aco3((fr )cos(^  d)  +  A  sin((}i)sin(^  d) 

=■  X  d  ♦  Y  cos(|^  d]  +  Z  sin(|^  d] 

where  D  =  milligal  equivalent 

d  =  counter  reading  (observable) 

Y  =»  Acos(j) 

adjusted  parameters 

Z  -  Asinij) 


The  condition  equation  is: 
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(31n(|l  dj  -  31n(f  d^^^))  2  .  0 


...  _('2Tr 


or 


ACY+ASZ+W-0  (W=xAd-W) 

o 


AC  -  [co3(f^  ]  -  C03(|i  )] 

j  j 


'2ir 


AS  »  [3in(|^  d^^  ]  -  3in(f^  d^^j.  )] 

j  J 


..._(•  2ir 


AC,  AS 


1  ““1 


AC.,  AS, 


AS„ 
r  r 


B  »  X 


1 

0 


-1 


-1 

0 


-1 


2r 


BB^  »  2x^1  (BB^)  ^  »  ~r  I  ;  (  I  -  unit  matrix) 

r  r  ;  2x^  ^  ^ 


T  T  -1  1 

N  -  A  (BB  )  A  -  — : 


2x 


ZAC  Z(ACAS) 

.2 


Z(ACAS)  ZAS 


-1 


2x 


ZAC^  ZAS^  -  (ZAS  aC)^ 


ZAS  -Z(aCAS) 


-Z(ACAS)  ZAC 


-  2a  x^  [ZAS^  ♦  ZAC^j 

TR(Z^)  -  TR(N  )  -  —  , - rTTm 


ZAS  j;aC  -  (EAS  ZAC) 
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•  •  '  »  * 


Let  0  -  d  „  ;  0  -  d  ;  0^-1  (0,  +  0,) 

IT  off  2Ton  tn21  2 


A0  -  2  ^®2  “  ®i  '  “  GOi^tant 


(Because  d  -  d  »  W/x) 
on  off 


2  2  2  2  2 
EAS  -  E(3in0_  -  3in0, )  -  E(2  0030  3inA0)  =  ^  3in  A0  E  cos  0 

2  1  m  m 


2  2  2  2  2 
EAC  -  E(co30_  -  0030.)  =  I(-23in0  3inA0)  »  4  sin  A0  E  sin  0 

2  1  tn  m 


r  r  r 

2  2  2  2  2  2 
E  AS  +  E  AC  -  4  sin  A0  E  (sin  0  +  cos  0  )  =  4  sin  A0  n 

11  ,  m  m 


E  AS^  EAC^  •>  (4  sin^A0  E  CO3^0  )  (4  3in^A0  E  3in^0  ) 

m  m 


4  2  2 

16  sin  A0  (E  cos  0  E  sin  0  ) 
m  m 


16  3in^A0  iE[^  (1  +  00320^)]  e[^  (1  -  CO320^)]} 


4  3in  A0  E  (1  +  cos20  )  E  (1  -  cos20  ) 
1  "  1 


4  sin  A0  (r  +  E  co320  )  (r  -  E  cos20  ) 

m  m 


4  2  2 

4  sin  A0  [r  -  (E  CO320  )  ] 

tn 
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2  2  2 
(IAS  AC)  =  (I  4  sin  A0  sin0_,  CO30^) 

m  m 


4  /  ^  1  \‘i 

1 6  sin  A0  ( I  -r  3in20  j 
2  in 


4  ,  2 

4  sin  A0  [l  3in20„) 

1  ” 


Hence:  IAS^  IAC^  -  (IAS  AC)^  =  4  3in**A0  [r^  -  (I  cos20  -  (I  3in20  )^] 

m  m 


2  2 
IAS  +  lAC 


2 

4  sin  A0  r 


lAS^  lAC^  -  (IAS  AC)^  4  sin^A0  Cr^  -  (I  00326  -  (I  3in20  )^] 

m  m 


3in^A0  [r^  -  (I  00320^)^  -  (I  3in20^)^] 


(note:  I  denotes  I  ) 


when  taking  6A  and  64)  as  adjusted  parameters,  similarly,  the  formula  for 

TR(Ijj), 


Tr(I 


0  ^  x^[(1  +  A^)r  *  (A^  -  1)  I  COS20  ] 
o  .  ,  m . 

)  _ ill _ i _ 

A‘^sin^A0  {r'^  -  [(  I  °°320  ]  +  (  I  sin20  )  ]} 

J  =  1  J  J  =  1  j 


can  be  also  derived  by  the  above  procedure. 
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APPENDIX  C 


Observations  and  d^^.^  (in  counter  units) 
for  the  four  numerical  examples  reported  in  Section  4.3.1 


Example 

1 

Example 

2 

‘^off 

^on 

‘^off 

^on 

1 . 

499.9969 

671 .6531 

499.9969 

671  .6531 

2. 

679.9990 

851.6742 

671 .6490 

843.3302 

3. 

859.9995 

1031 .6756 

843.2995 

101 4.9977 

4. 

1039.9965 

121 1.7094 

1014.9465 

1186.6433 

5. 

1220.001  4 

1391 .6961 

1 186.601 4 

1358.31 38 

6. 

1400.001  4 

1571  .6943  * 

1 358.2514 

1529.9426 

7. 

1579.9985 

1751 .6857 

1529.8985 

1701 .5983 

8. 

1759.9991 

1931 .6956 

1701 .5491 

1873.2659 

9. 

1939.9988 

211  1  .6868 

1873.1988 

2044.8900 

10. 

2120.0022 

2291 .7102 

2044.8522 

2216.5620 

1 1  . 

2300.0013 

2471 .6957 

2216.5013 

2388.1869 

12. 

2479.9995 

2651.6749 

2388.1495 

2559.8426 

13. 

2660.0017 

2831  .6626 

2559.8017 

2731 .4715 

14. 

2839.9986 

3011.6432 

2731 .4486 

2903.1125 

15. 

3019.9993 

3191  .6340 

2903.0993 

3074.7430 

16. 

3200.0012 

3371 .6290 

3074.7512 

3246.3835 

17. 

3379.9988 

3551 ,6326 

3246.3988 

3418.0432 

18. 

3559.9972 

3731.5970 

3418.0472 

3589.6530 

19. 

3739.9999 

391 1 .6010 

3589.6999 

3761  .3295 

20. 

3920.0012 

4091 .5772 

3761 .3512 

3932.9417 

21 . 

4100.0013 

4271 .5766 

3933.0013 

4104.5947 

22. 

4279.9995 

4451.581 9 

4104.6495 

4276.2454 

23. 

4460.0024 

4631 .6072 

4276.3024 

4447,8740 

24. 

4640.0018 

481 1.5961 

4447.9518 

4619,5475 

25. 

4819.9974 

4991  .6133 

4619.5974 

4791  .2059 

26. 

5000.001 1 

5171 .61 40 

4791 .2511 

4962.8673 

27. 

5179.9974 

5351 .6083 

4962.8974 

5134.5119 

28. 

5360.0023 

5531 .6340 

5134.5523 

5306. 1747 

29. 

5539.9978 

571 1 .6491 

5306.1978 

5477.81 72 

30. 

5719.9976 

5891  .6616 

5477.8476 

5649.4904 
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